THE QUANTUM CONTENT OF THE GLUING EQUATIONS 
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Abstract. The gluing equations of a cusped hyperbolie 3-nianifold M are a system of 
polynomial equations in the shapes of an ideal triangulation T of M that describe the 
complete hyperbolic structure of M and its deformations. Given a Neumann-Zagier datum 
(comprising the shapes together with the gluing equations in a particular canonical form) we 
define a formal power series with coefficients in the invariant trace field of M that should (a) 
agree with the asymptotic expansion of the Kashaev invariant to all orders, and (b) contain 
the nonabelian Reidemeister-Ray-Singer torsion of M as its first subleading "1-loop" term. 
As a case study, we prove topological invariance of the 1-loop part of the constructed series 
and extend it into a formal power series of rational functions on the PSL(2, C) character 
variety of M . We provide a computer implementation of the first three terms of the series 
using the standard SnapPy toolbox and check numerically the agreement of our torsion 
with the Reidemeister-Ray-Singer for all 59924 hyperbolic knots with at most 14 crossings. 
Finally, we explain how the definition of our series follows from the quantization of 3d 
hyperbolic geometry, using principles of Topological Quantum Field Theory. 
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1. Introduction 

1.1. The Kashaev invariant and perturbative Chern-Simons theory. The Kashaev 
invariant (K) at G C of a knot K in 3-space (for A^ = 2,3,...)isa powerful sequence of com- 
plex numbers determined by the Jones polynomial of the knot [Jon87] and its cablings [Tur88, 
Wit89]. The Volume Conjecture of Kashaev and Murakami- Murakami [Kas97, MMOl] re- 
lates the Kashaev invariant of a hyperbolic knot K with the hyperbolic volume Vol(M) of 
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its complement M = \ K [Thu77] 

A generalization of the Volume Conjecture [Guk05] predicts a full asymptotic expansion of 
the Kashaev invariant to all orders in 1/iV 

(1.2) {K)M''^'^ZM{2in/N) 
for a suitable formal power series 

(1 3 \ 2TTi 

n>2 ^ 

The formal power series Zm{K) in (1-3) is conjectured to coincide with the perturbative par- 
tition function of Chern-Simons theory with complex gauge group SL(2, C) along the discrete 
faithful representation po of the hyperbolic manifold M. Combining such an interpretation 
with further conjectures of [DGLZ09, Gar08, GLll] one predicts that 

• Sm,o = i(yolM + ^CSa/) G C/(47r^Z) is the complexified volume of M {cf. [Thu82, 
Neu92]). 

• Sm,i should be related [Wit91, BNW91, GM08] to the nonabelian Ray-Singer torsion 
[RS71], which ought to equal (c/. [Miil93]) the combinatorial nonabelian Reidemeister 
torsion. More precisely, [DC, Conj.1.8] we should have 

(1.4) Tf,=47T'exp{-2SM,i)eEl, 

where r]^ is the nonabelian Reidemeister-Ray-Singer torsion of M with respect to 
the meridian [Por97, Dub06], and Em is the invariant trace field of M. 

• For n > 2, the n-loop invariants SM,n are conjectured to lie in the invariant trace 
field Em [DGLZ09, Gar08]. 

The generalization (1.2) of the Volume Conjecture has been numerically verified for a few 
knots using either state-integral formulas for Chern-Simons theory when available [DGLZ09] 
or a numerical computation of the Kashaev invariant and its numerical asymptotics, lifted 
to algebraic numbers [GZ, Gar]. 

Our goal is to provide an exact, combinatorial definition of the formal power series Zm{K) 
via formal Gaussian integration using the shape parameters and the Neumann- Zagier ma- 
trices of a regular ideal triangulation of M. Our definitions 

• express the putative torsion exp(— 25'm,i) and the n-loop invariants SM,n manifestly 
in terms of the shape parameters Zi and the gluing matrices of a regular ideal trian- 
gulation T of M; 

• manifestly deduce that the putative torsion and the n-loop invariants for n > 2 are 
elements of the invariant trace field; 

• explain the difference of Zm{K) for pairs of geometrically-similar knots studied by 
Zagier and the second author; 

• provide an effective way to compute the n-loop invariants using standard commands 
of the SnapPy toolbox [CDW] — as demonstrated for n < 3 for hyperbolic knots with 
at most 14 crossings; and 
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• allow efficient tests of the asymptotics of the Volume Conjecture (1.2), the "1-loop 
Conjecture" (1.4) and other conjectures in Quantum Topology. 

We note that we only define exp(— 2S'm,i) up to a sign, and Sm,2 modulo Z/24. All higher 
n-loop invariants are defined unambiguously. 

Although we give a purely combinatorial definition of Zm{K) without any knowledge of 
state integrals or Chern-Simons theory with complex gauge group, in Section 5 we explain 
how our definition follows from the state-integral model of the ffist author [Dim] and its 
perturbative expansion. 

1.2. The Neumann-Zagier datum. All manifolds and all ideal triangulations in this pa- 
per will be oriented. The volume of a hyperbolic manifold M, appearing in the Volume 
Conjecture and contributing to Sma-, is already known to have a simple expression in terms 
of shape parameters of a regular ideal triangulation, i.e., one that recovers the complete 
hyperbolic structure of M. (For extended discussion on regular triangulations, see Section 
4.) If T = { Aj}^^ is a regular ideal triangulation of M with shape parameters Zi G C\ {0, 1} 
for i = 1, . . . , iV, then (c/. [DS82, NZ85]) 

N 

(1.5) Vol(M) = ^D(z,), 

1=1 

where D[z) := Im(Li2(2)) +arg(l — z) log \z\ is the Bloch-Wigner dilogarithm function. This 
formula can also be interpreted as calculating the image of the class [M] := XIJ-^*] £ of 
M in the Bloch group B under the natural map £) : ,B — t- M. An analogous formula, using 
the class of M in the "extended" Bloch group gives the full complexified volume Sma 
[Neu92, Neu04, GZ07, Zic09]. 

It is natural to ask whether the class of M in S determines not only Sma but the higher 
SM,n as well. This question was posed to the authors several years ago by D. Zagier. Sub- 
sequent computations [GZ, Gar] indicated that a positive answer was not possible. For 
example, there is a family of pairs of pretzel knots ((-2,3,3 + 2p), (-2,3,3 - 2p)) for 
p = 2, 3, as well as the figure-eight knot and its sister, which all have the same class in the 
Bloch group (and classes differing by 6-torsion in the extended Bloch group), but different 
invariants SM,n for n > 1. 

The extra information necessary to determine the Sm^u can be described as follows. Recall 
that if T is a regular ideal triangulation of M with N tetrahedra, its shapes z = (2:1, ... , zn) 
satisfy a system of polynomial equations, one equation for every edge, and one imposing 
parabolic holonomy around the meridian of the cusp [Thu77, NZ85]. Let us set 

(1.6) z[ = {l- z,)-' , = 1 - Z-' . 
The equations can then be written in the form 

N N 

(1.7) n ^/"^^^'^ = n (1 - ^j'T^' = ±1 ' = 1' ^) ' 

i=i j=i 

where A and B are N x N square matrices with integer entries, which we call the Neumann- 
Zagier matrices following [NZ85]. 
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Definition 1.1. If T is a regular ideal triangulation of M, its Neumann- Zagier datum (resp., 
enhanced Neumann- Zagier datum) is given by the triple 

f3r={z,A,B) resp., /3r = A, B, /) , 

where / is a flattening of T, a collection of integers that we define in Section 2. 

As we will discuss in detail in Section 2, implicit in the above definition is the dependence 
of f3j- and Pj- on the following choices: 

(1) a pair of opposite edges for every oriented ideal tetrahedron (a so-called choice of 
quad type), 

(2) an edge of T, 

(3) a meridian loop in the boundary of M in general position with respect to T, 

(4) a flattening. 

1.3. The 1-loop invariant. 

Definition 1.2. Given a one-cusped hyperbolic manifold M with regular ideal triangulation 
T and enhanced Neumann- Zagier datum Pj- we define: 

(1.8) rr := ±^ det (AA,,, + BA~^)z^"z"~^ e Em/{±1} , 

where A^ := diag(2i, z^) and A^// := diag(2;^', z'J^) are diagonal matrices, and z^" z"~f : = 

Note that tj- takes value in the invariant trace field of M and is only defined up to a sign. 
In Section 3 we will show that 

Theorem 1.3. Tq- is independent of the quad type ofT, the chosen edge ofT, the choice of 
a meridian loop, and the choice of flattening. 

We now consider the dependence of r-j- on the choice of a regular ideal triangulation of M. 
It is well known that the set X of ideal triangulations of a cusped hyperbolic manifold is 
non-empty [Cas65] and connected by 2-3 moves [Mat87, Mat07, PieSS]. That is, a sequence 
of 2-3 moves can be used to take any one ideal triangulation to any other. The subset Xp^ 
of X of regular triangulations is also non-empty, see Section 4. Topologically, these are the 
triangulations without any univalent edges [Cha03, BDRV, DGll, Til]. We will prove in 
Section 3 that 

Theorem 1.4. t-j- is constant on every component of Xp^ connected by 2-3 moves. 

1.4. Expectations. We may pose some questions and conjectures about the 1-loop invariant 
rr and the structure of the set Xp^y Let us begin with two questions whose answer is 
unfortunately unknown. 

Question 1.5. Is Xp^^ connected by 2-3 moves? 
Question 1.6. Is tj- constant on the set Xp^^l 
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Clearly, a positive answer to the first question implies a positive answer to the second. 

Despite unknown answers to the above questions, with additional effort we can still define 
a distinguished component of A'pg , and thus obtain a topological invariant of M. Namely, let 
Xep C A'pg denote the subset that consists of regular refinements of the canonical (Epstein- 
Penner) ideal cell decomposition of M [EP88]. Xep is a finite set and generically consists of 
a single element. In Section 4.2 we will show the following. 

Proposition 1.7. Xep lies in a connected component of Xp^. Consequently, the value of t-j- 
on Xep is a topological invariant tm of M. 

Admittedly, it would be more natural to show that rr is constant on all of Xp^^ . Proposition 

1.7 appears to be an artificial way to construct a much needed topological invariant of cusped 
hyperbolic 3-manifolds. 

Our next conjecture compares our torsion tm with the nonabelian Reidemeister torsion 
of M with respect to the meridian defined in [Por97, Dub06]. 

Conjecture 1.8. For all hyperbolic knot complements we have = ±rA/. 

Numerical evidence for the above conjecture is presented in Appendix D using Dunfield's 
computation of rfj via SnapPy [DP J]. Observe that both sides of the equation in Conjecture 

1.8 are algebraic numbers (defined up to a sign) that are elements of the invariant trace 
field of M. Moreover, if M has a regular ideal triangulation with N ideal tetrahedra and its 
fundamental group is generated with r elements, then tm and rj[| are essentially given by 
the determinant of square matrices of size and 3r — 3, respectively. It is still unclear to 
us how to relate these two matrices or their determinants. 

By definition, rj[^ G E*j^. Thus, a mild but important corollary of Conjecture 1.8 is that 
Tm is nonzero. This is a crucial ingredient, necessary for the definition of the higher loop 
invariants SM,n using perturbation theory. 

1.5. The higher-loop invariants. In this section we define the higher loop invariants Sq-^n 
for n>2. They are analyzed in detail in Section 5, using a state integral (5.2). The result, 
however, may be summarized as follows. Let us introduce a formal power series 



where Bn is the n^^ Bernoulli number (with Bi = +1/2), and Lm{z) is the m**^ polylogarithm. 
Note that Lim(2) G {1 — z)~^~^'Z[z] is a rational function for all nonpositive integers m. This 
formal series arises from the asymptotic expansion of the quantum dilogarithm function 
[Bar99, FK94, Fad95], which is the Chern-Simons partition function of a single tetrahedron. 

We fix an enhanced Neumann- Zagier datum Pj- = [z, A, B, /) of an oriented 1-cusped 
manifold M and a regular ideal triangulation T with A^ tetrahedra. Let u = Af + B/". We 
assume that 



The condition det(B) ^ can always be satisfied with a suitable labeling of shapes (Lemma 
A. 3). It follows (c/. Equation (2.9) below) that 



(1.9) 




det(B) 7^ 



(1.10) 



n 
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with A^/ = diag(2;^, ...,2;^), is a non-degenerate symmetric matrix. We define 

(1.11) fT,ni^;z) = exp ( - —x^R-V + -/^B-^A/ ] n^r.(x„z,) G Q(^)[x,;i5] , 

^ ^ j=i 

where Y and z = (^1, z^). A formal power series fhix) G Q(2;)|x, /iz] has 

a. formal Gaussian integration, given by (c/. [BIZ80]) 

J axe -2 

This integration is defined by expanding fh{x) as a series in x, and then formally integrating 
each monomial, using the quadratic form "H"^ to contract x-indices pairwise. 

Definition 1.9. With the above conventions, we define 

(1.13) exp 5r,„(^)r"i j := {fT,H{x- z)) . 

Notice that the result involves only integral powers of h. Moreover, 57-^^ G r^'^"'^'^Q[z, z', z"\ 
for all n > 2. Thus, we can also write 

n=2 / n=\ 7" 

where S-j^n G Qf-z, 2', -z"]. 

1.6. Feynman diagrams. A convenient way to organize the above definition is via Feynman 
diagrams, using Wick's theorem to express each term S-y^n as a finite sum of connected 
diagrams with at most n loops. This is well known and explained in detail, e.g. in [HKK"'"03, 
Ch.9], [BIZ80, P0IO5]. 

The Feynman rules for computing the Sr,n-, described in Section 5, turn out to be the 
following.^ One draws all connected graphs D with vertices of all valencies, such that 

(1.15) -^(-D) := (# 1-vertices + # 2-vertices + # loops) < n. 
In each diagram, the edges represent an A x A propagator 

(1.16) propagator : 11 = hH"^ , 

(k) 

while each /c- vertex comes with an A- vector of factors F- , 

P = ak I ~ 

where Ofc = 1 (resp., 0) if = 1,2 (resp., k > 3). The diagram D is then evaluated by 

(k) 

contracting the vertex factors T- with propagators, multiplying by a standard symmetry 



^To derive these from (1.12), one should first rescale x ^ h 2; 
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factor, and taking the part of the answer. In the end, SM,n is the sum of evaluated 
diagrams, plus an additional vacuum contribution 

TV 




ra > 3 

To illustrate the above algorithm, we give the explicit formulas for 5*2 and S3 below. 

1.7. The 2-loop invariant. The six diagrams that contribute to Sm,2 are shown in Figure 
1, together with their symmetry factors. 




-<3 U 



Figure 1. Diagrams contributing to S'a/,2 with symmetry factors. The top 
row of diagrams have exactly two loops, while the bottom row have fewer loops 
and additional 1-vertices and 2-vertices. 



Their evaluation gives the following formula for S'7-,2: 



:i.i9) 



S- 



r,2 



coeff 



-F^^Vn- )^ + -n- F';"'^!!- F^.'^^n - + — F'^'^^fn 



.(3) 



(3)i 



^(3), 



12 



\3p(3) 



J F-^^n •F^-''TI 

^2 J 



^(3)i 



+ fw 



where all the indices i and j are implicitly summed from 1 to and coeS[f {h) , h] denotes the 
coefficient of of a power series f{h). Concretely, the 2-loop contribution from the vacuum 
energy is 

F(o) = l/B-^A/ - ^'v and the four vertices that appear only contribute at 

leading order, 
(1.20) 



(B-V), 



.(2) 



Ziz: 



12 



.(4) 



Zi{l + Zi)z: 



/3 



2 ' * 2 ' * h ' ' h 

We expect 5'7-,2 to be well defined modulo Z/24, and this is exactly what happens in hundreds 
of examples that we computed. 

1.8. The 3-Ioop invariant. For the next invariant 5*7-^3, all the diagrams of Figure 1 con- 
tribute, collecting the coefficient of fi^ of their evaluation. In addition, there are 34 new 
diagrams that satisfy the inequality (1.15); they are shown in Figures 2 and 3. Calculations 
indicate that the 3-loop invariant 6*7-^3 is well defined, independent of the regular triangula- 
tion T. The invariants S-Yfl-, Tr, Sffi have been programmed in Mathematica and take 
as input a Neumann- Zagier datum available from SnapPy [CDW]. 
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For the 4- loop invariant, there are 291 new diagrams. A python implementation will be 
provided in the future. For large ra, one expects about nl^C"' diagrams to contribute to Sn- 




Figure 2. Diagrams with three loops contributing to 



o^— o shI) i-Q 



;<XD 




SO 




1 
4 



Figure 3. Diagrams with 1-vertices and 2- vertices contributing to S3. 



1.9. Generalizations. There are several natural extensions of the results presented above. 
First, one could attempt to prove the independence of the all-loop invariants Zf{h), including 
the entire series of S^-^nS, under 2-3 moves and different choices of Neumann-Zagier datum. 
This was done non-rigorously in [Dim], but a full mathematical argument in the spirit of 
Theorems 1.3 and 1.4 is still missing. We hope to address this in future work. 
In a different direction, one can extend the formulas for tj- and S']-,n to 

- manifolds with multiple cusps, 

- representations other than the discrete faithful. 
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- representations with non-parabolic meridian holonomy, 

- non-hyperbolic manifolds. 

The only truly necessary condition is that a 3-manifold M have a topological ideal triangu- 
lation T that — upon solving gluing equations and using a developing map — reproduces 
some desired representation p : tti{M) — )■ PSL(2,C). We call such an ideal triangulation 
p-regular, and in Section 4 we will briefly discuss most of the above generalizations. In 
particular, we will demonstrate in Sections 4.6 and 5.5 how to extend t-j-, S-j-^n to rational 
functions on the character variety of a (topologically) cusped manifold. The generalization to 
multiple cusps is also quite straightforward, but left out mainly for simplicity of exposition. 

Acknowledgement. The authors wish to thank Nathan Dunfield, Walter Neumann, Jose- 
phine Yu and Don Zagier for many extremely enlightening conversations. 

2. Mechanics of triangulations 

We begin by reviewing the gluing rules for ideal hyperbolic tetrahedra and the equations 
that determine their shape parameters. We essentially follow the classic [Thu77, NZ85], but 
find it helpful to work with additive logarithmic (rather than multiplicative) forms of the 
gluing equations. Recall that all manifolds and all ideal triangulations are oriented. 

2.1. Ideal tetrahedra. Combinatorially, an oriented ideal tetrahedron A is a topological 
ideal tetrahedron with three complex shape parameters {z, z', z") assigned to pairs of opposite 
edges (Figure 4). The shapes always appear in the same cyclic order (determined by the 
orientation) around every vertex, and they satisfy 



(2.1a) zz'z" = -l, 

and 

(2.1b) z" + z-^-l = Q. 

In other words, z' = — z) and z" = 1 — z~^. We call the tetrahedron non- degenerate if 
none of the shapes take values in {0, 1, oo}, i.e., z, z\ z" G C*\{1}. It is sufficient to impose 
this on a single one of the shapes. 



(oo) 




Figure 4. An ideal tetrahedron 
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Borrowing common terminology from the theory of normal surfaces, cf. [Bur, Kan05, 
KR04, Til08], we define the quadrilateral type (in short, quad type) of A to be the distin- 
guished pair of opposite edges labelled by z. Clearly, there is a three-fold choice of quad type 
for any oriented ideal tetrahedron. Different choices correspond to a cyclic permutation of 
the vector {z, z,' , z"), which leaves relations (2.1) invariant. 

Geometrically, the shape parameters determine a PSL(2, C) structure on A. Equivalently, 
they determine a hyperbolic structure, possibly of negative volume. We can then describe 
the ideal hyperbolic tetrahedron A as the convex hull of four ideal points in hyperbolic three- 
space H^, whose cross-ratio is z (or z', or z"). Each shape z fixes the complexified dihedral 
angle on the edge it labels, via 

(2.2) z = exp (torsion + i angle) , 

and similarly for z', z" . We will choose logarithms (Z, Z', Z") such that (2;, z\ z") = (e^, e^', e^") 
and 

(2.3) Z + Z' + Z" = m . 

This condition has the natural interpretation that the sum of dihedral angles going around 
any vertex is vr, and that the torsion vanishes. 

Recall that the hyperbolic volume of an oriented ideal tetrahedron with shape z is given 
by the Bloch-Wigner function D{z), cf. (1.5) [DS82, NZ85]. This function satisfies D{z) = 
D{z') = D{z"), consistent with the cyclic permutation symmetry of A. The tetrahedron has 
positive volume if and only if '^{z) > 0. In this case, taking standard logarithms {Z, Z', Z") 
of the shapes (such that < '^{Z) < vr, etc.) automatically implies (2.3). 

2.2. Gluing. Let M be an oriented one-cusped manifold with an ideal triangulation T = 
{Aj}^;^ and a choice of quad type. Assign shapes {zi, z[, z'l) (and logarithms (Zj, Z[, Z'D) to 
each tetrahedron Aj according to the choice of quad type and cyclic order induced by the 
orientation of M. A simple Euler characteristic argument shows that the triangulation has 
exactly N edges Ej, I = 1, N. Each edge Ej gives rise to a gluing equation of the form 

N 

(2.4) Ej: ^ (G,,Z, + G',,Z: + G';,Zf ) = 27rz I=1,...,N, 

1=1 

where G/j G {0, 1, 2} (resp., G'^^, G"J is the number of times an edge of tetrahedron Aj with 
parameter Zi (resp., z', z'l) is incident to the edge Ei in the triangulation. 

The Equations (2.4) are not all independent. For a one-cusped manifold, every edge begins 
and ends at the cusp, which implies Yl!i=i — Tld=i ^'u — S/Li ^'u = 2, and therefore 
that the sum of the left-hand sides of Equations (2.4) equals 2'iriN. This is the only linear 
dependence. In general, there is one relation per cusp of M. 

When the hyperbolic structure of M is complete, the PSL(2, C) holonomy of the meridian^ 
fi is parabolic, ~ ( { )• This imposes an additional linear constraint on the shapes, of the 



If M is not a knot complement, a canonical peripheral pair of a meridian and longitude may not exist. 
A meridian cycle must then be chosen, and held fixed. The torsion and invariants 5j\/_„ will depend on the 
choice of a peripheral system in a predictable way, cf. [Yam08] . 
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form 

N 

(2.5) /i : ^ (GAT+i^jZj + G'^^i jZ- + G'l^^i iZ'/) = , 

i=l 

with Gat+i^j, G^_,_i j, G'l^_^_i j G Z. To calculate these integers, one can draw the 2-dimensional 
torus boundary of the cusp neighborhood, triangulated by the intersections of horospheres 
with tetrahedron vertex neighborhoods (c/. [NZ85] and Section 2.5 below). Then one should 
choose a normal path in this torus that is homotopic to the meridian and has no self- 
intersections. The equation (2.5) is obtained by adding and subtracting the logarithmic 
shape parameters on the angles subtended by the path, with signs corresponding to orienta- 
tion [Neu92]. The convention for choosing vs. '— ' signs is immaterial in (2.5), since the 
right-hand side vanishes. It turns out that the logarithmic constraint (2.5) is independent 
of the precise choice of meridian loop — modulo the logarithmic edge constraints (2.4) and 
the relation (2.3). 

2.3. The Neumann- Zagier matrices. The linear conditions (2.3) (2.4), and (2.5) have 
both symmetry and redundancy. We have already observed that any one of the edge con- 
straints (2.4) can be removed. Let us then ignore the edge I = N. We can also use (2.3) 
to eliminate one of the three shapes for each tetrahedron. We choose this canonically to be 
though which pair of edges is labelled Z'^ depends on the choice of quad type for the 
tetrahedron. Then the first — 1 edge equations and the meridian equation are equivalent 
to 

N 

(2.6) Yl (^^^-^i + ^i^^j) =^^^^^ i = l,...,N, 
where 

(2.7) = 1 ' . Bi, = l 

I Gat+ij - G^_^ij i - N I G^_^ij - G'^_^i_j- i - N; 



and 
(2.8) 



We will generally assume Z, Z", and u to be column vectors, and write A Z + B Z" = mv . 
The data (A, B, u) (in fact, the full logarithmic equations (2.4)-(2.5)) can be obtained from 
SnapPy [CDW], illustrated in in Appendix D. 

The Neumann- Zagier matrices A and B have a remarkable property: they are the top two 
blocks of a 2N x 2N symplectic matrix [NZ85]. It follows that 

(2.9) AB^ = BA^, 

and that the x 2N block (A B) has full rank. This symplectic property is crucial in 
defining the state integral of [Dim], in defining our formal power series invariant Zm{K)i and 
in the combinatorial proofs of topological invariance of the 1-loop invariant. 
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If we exponentiate the equations (2.6) and combine them with the nonhnear relation (2.1b) 
for each tetrahedron, we obtain 

(2.10) z^^"^ = ^^(1 - z-^)^ = (-1)% 

where := W- z'^'^ . These equations in variables fully capture the constraints imposed 
by the gluing. For hyperbolic M, a triangulation T is regular precisely when one of the 
solutions to (2.10) corresponds to the complete hyperbolic structure. This "geometric" 
solution is distinguished as the only solution with all '^{zi) > 0, i.e., with all tetrahedra 
having non-negative volume. 

2.4. Flattening. Given an ideal triangulation T of M, a flattening [Neu92, Neu04] is a 
collection of 3N integers (/j, //, //') that satisfy 

(2.11a) /^ + /; + /;' = i, z = i,...,iv, 

N 

(2.11b) (G/^/^ + G/J/ + G/JD = 

i=l 

Equivalently, it is sufficient to impose (2.11a) and 

(2.12) Af + Bf" = u. 

Evidently, a flattening is a discretized version of logarithmic shape parameters, satisfying 
only the linear gluing equations. It is so named because it provides a way to "flatten" 
tetrahedra: if {Zi, Z") satisfy (2.3) and (2.4), then the flattened shapes 

(Z„ Zi Z':) = {Z,- ^7^/^, Z[ - mfi Z'l - n^fl) 

satisfy the same equations with the right-hand sides set to zero. 

Remark 2.1. By a theorem of [Neu92], a flattening always exists. 

We note that given an enhanced Neumann-Zagier datum (2, A, B, /), we can easily recon- 
struct the integer vector v via (2.12). 

2.5. Example: 4i. As an example, we describe the enhanced Neumann-Zagier datum of 
the figure-eight knot complement M. It has a well known regular ideal triangulation T 
consisting of A = 2 tetrahedra, to which we assign logarithmic shape parameters (Z, Z', Z") 
and {y^,W\W"). 




2 / = 1,...,A 

/ = A + 1 



Figure 5. The boundary of the cusp neighborhood for the figure-eight knot. 
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A map of the boundary of the cusp neighborhood is shown in Figure 5. We have chosen 
one of 3^ possible cychc labehngs by Z's and W's. Each of the edges intersect the cusp twice, 
so it is easy to read off from Figure 5 that the edge constraints (2.4) are 

El : 2Z + Z" + 2W + W" = 2TTi 

E2 : 2Z' + Z" + 2W' + W" = 2TTi . 

The sum of the left-hand sides is automatically 47ri, so we can choose to ignore the second 
constraint. If we choose the meridian path fi as in Figure 5, the meridian constraint (2.5) is 

/i : -Z' + W = 

Putting together the first edge constraint and the meridian into matrices, we have 
2 2\ f Z\ fO 0\ f Z'\ fl l\ f Z'\ . (2 



\) \ W ) ' \-\ OJ \W'J \o oj \W'J \0 
and after using Z + Z' + Z" = W + W + W" = in to eliminate Z' and W, we get 

2 2\ f Z\ fl l\ f Z'\ . (2 



1 \) \ W ) ' \\ Oj \W'J \1 
From this last expression, we can read off 

(2.13) ^=(ll)' ''^Go)' "=(1 

The two gluing equations (2.10) are then 

(2.14) z'^w'^z"w" = 1 , zwz" = -1 , 

with z" = 1 — and w" = \ — w~^. The solution for the complete hyperbolic structure is 
z = w = e™l'^ . 

Finally, a flattening (/^, /^'; /^) e i.e. an integer solution to A/ + B/" = v 

and / + /' + /" = 1, is given by 

(2.15) (/„ fl U, C a = (0, 1, 0; 1, 0, 0) 

It is easy to see that every flattening has the form {a,b,l — a — b;b,a,l — a — b) for integers 
a, b. 

3. Topological invariance of our torsion 

Given a one-cusped hyperbolic manifold M with regular triangulation T = {Aj}^^ and 
Neumann- Zagier datum /3j- = {z, A, B, /), we have proposed that the nonabelian torsion is 
given by 

(3.1) rr := ±^ det (AA,. + BA;^)z^"z"-^ e Em/{±1} , 

where A^ = diag(2;i, ^at), and similarly for A^". Since {z,z',z") G Em we must have 
Tj- G Em as well. 

After a brief example of how the formula (3.1) works, we will proceed to prove Theorems 
1.3 and 1.4 on the topological invariance of Tf. We saw in Section 2 that the Neumann- Zagier 
datum depends not only on a triangulation T, but also on a choice of 
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(1) quad type for T, 

(2) one edge of T whose gluing equation is redundant, 

(3) normal meridian path, 

(4) flattening /. 

We will begin by showing that tj- is independent of these four choices, and then show that 
it is invariant under 2-3 moves, so long as the 2-3 moves connect two regular triangulations. 

The four choices here are independent, and can be studied in any order. However, in 
order to prove independence of flattening, it is convenient to use a quad type for which the 
matrix B is non-degenerate. Such a quad type can always be found (Lemma A. 3), but is 
not automatic. Therefore, we will first show invariance under change of quad type, and then 
proceed to the other choices. It is interesting to note that of all the arguments that follow 
(including the 2-3 move), independence of flattening is the only one that requires the use of 
the full gluing equations z^z"^ = (—1)''. 

3.1. Example: 4i continued. To illustrate the Definition 1.8, consider the figure-eight 
knot complement again. From Section 2.5, we already have one possible choice for the 
Neumann- Zagier matrices (2.13) and a generic flattening (2.15). We use them to obtain 



idet 



2 2 
1 1 



2w" 
1 

2^" 



■det 



■det 



2z" + z-^ 
z" + z-^ 

z" + 1 





w" 
2w" 



w 



1 1 
1 

1 











w 



w 



."-1 



w 



w 



w 



(3.2) 



- 2 
2^ 



w 



-3, 



1 



w ' 



w 



-1 



where at intermediate steps we used z" + z~ 
substituted the discrete faithful solution z = z" = w = w" = e 

The invariant belongs to the invariant trace field i?4^ 
torsion of the figure-eight knot complement [DFJ]. 



in/ 3 



= 0, and at the end we 
-3), and agrees with the 



3.2. Independence of quad type. Now, let us fix a manifold M, a triangulation T with 
tetrahedra, and an enhanced Neumann- Zagier datum Pj- = {z, A, B, /). 
To prove independence of quad type, it is sufficient to check that rr is invariant under a 

cyclic permutation of the first triple of shape parameters (^i, z[, z'{), while holding fixed the 



1^1, 



.Zn), U = (l/i. 



choice of meridian loop and redundant edge. Let us write z 
f = (/i,--,/7v)^, and 

(3.3) A = («!, 02, ■ ■ ■ ,aAr) , B = 62, ■ ■ ■ ,&7v) , 

in column notation. After the permutation, a new Neumann- Zagier datum is given by 
[z, A, B, /) where 

(3.4) z = {z[, Z2, zn) , z* = {z'l, Z2, z'jsf) , z" = {zi, Z2, z'^) , 



(3.5) A = (-61, 



^2, 



qn) , B = (ai - bi, 62, ■ ■ ■ ,bN) 



V 



[nx - 61,77,2, 
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The new shapes satisfy z z' 
by permuting 
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-l)*^. We also naturally obtain a new flattening (/, /', /") 



(3.6) 



f — ifi: f2, In) 5 /' — (/I'l •■■5 /at) ! /" — /2 ; /^') j 



this is an integer solution to A/ + B/" = u and / + /' + /" = !. 

The torsion tj- (1.8) consists of two parts, a determinant and a monomial correction. By 
making use of the relations zi + z[~^ — 1 = and ziz[z'( = —1, we find the determinant with 
the permuted Neumann- Zagier datum to be 



det (AAj// + BA- ^) = det ( - biZi + (ai - bi)z[ \ a24' + ^2-^2 



det [aiz[ — hi, a^z'^ + 1)2'^ 



"1 ' 



(3.7) 



By simply using 
transforms as 



1 



-zi det (aiz^' + \ 02^' + ^222'^ ' ' ' 5 «iv^7v + ^Af^^v^) 
-zi det (AA^// + BA;^) . 

— 1 and /i + /( + /{' = 1, we also see that the monomial correction 



P"z"-^ 



fi'z'J-h 



Zi^i Z'{ 



(3.8) 



zf"z"~H 



zf"z"-f 



z. 



The extra factors in the two parts of the torsion precisely cancel each other, leading in 
the end to 



(3.9) 



det (AAj,, + BAri)#"i"-/ = (-l)/i+i det (AA,» + BA;^) 



zf"z"-f 



This is just as desired, showing that the torsion is invariant up to a sign. 

3.3. Independence of the choice of edge. We fix M, T, = (2, A, B,/). In order to 
choose matrices A, B, we must ignore the redundant gluing equation corresponding to an 
edge of T. This was discussed in Section 2.3. Suppose, then, that we choose a different edge 
to ignore. For example, if we choose the (A^ — 1)^* rather than the A^*^ (and keep the same 
quad type and meridian path), then we obtain new Neumann- Zagier matrices A, B, which 
are related to the original ones as 



(3.10) 


A = P(A._ 


1,N) A 


7 


B = P^N- 




B, 


where 


















/ 1 





... 





°\ 









1 


... 








(3.11) 


P{N-1,N) '■ = 








... 1 












-1 - 


-1 


. . . _i 


-1 


-1 











... 





1 / 
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Similarly, eliminating the rather than the A^"^ edge constraint is implemented by mul- 
tiplying with a matrix P(i,n) whose row is filled with — I's. Any such matrix satisfies 
detP(/,7v) = -1. 

In the formula for rr, only the determinant part is affected by a change of edge. Then 
(3.12) det (AA^// + BAj^) = det (P(7,jv)(AAy/ + BA;^)) = - det {AA,„ + BA;^) , 
leading to invariance of t-j-, up to the usual sign. 

3.4. Independence of meridian path. Recall that an ideal triangulation on M induces a 
triangulation of its boundary torus dM. Consider two simple closed meridian loops in dM in 
general (normal) position with respect to the triangulation of dM. Recall that these paths 
are drawn on the triangulated 2-dimensional torus dM where faces of tetrahedra correspond 
to edges in the 2-dimensional triangulation, and edges of tetrahedra to vertices. In particular, 
for a one-cusped manifold M, every edge of the triangulation intersects a pair of vertices on 
the boundary dM. 




Figure 6. The fundamental move for changing a meridian path. Here, we 
deform through an edge Ej with gluing constraint Xj = Z + W + Ri + ... + 
Rj + Si + ... + Sf = 27rz. 

We can deform one of our meridian paths into the other by using repeated applications of 
the fundamental move shown in Figure 6 — locally pushing a section of the path across a 
vertex of dM. Thus, it suffices to assume that the two paths only differ by one such move. 
Suppose that we cross the P^ edge (by Section 3.3 we may assume that / ^ N), which has 
a combinatorial gluing constraint 

N 

(3.13) Xj := ^ [GnZ, + G'^X + G'/.Zf) = 2m , 

i=l 

and that the two tetrahedra where the paths enter and exit the vicinity of the edge have 
parameters {Z, Z', Z") and (W, W, W"), as in the figure. We do not exclude the possibility 
that {Z,Z',Z") and (W,W' ,W") both coincide with the same triple (Zj, Z^', Zf), in some 
cyclic permutation. Then the difference in the logarithmic meridian equations (2.5) for the 
two paths will be 

(3.14) 6^ = ±{Xj - {Z + Z' + Z") -{W + W' + W")) . 
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Note that two logarithmic meridian constraints that differ by (3.14) are compatible and 
equivalent, since upon using the additional equations Xj = 27ri and Z + Z' + Z" = W + 
W + W" = iTT, we find that 5^ = 0. A discretized version of this observation demonstrates 
that that the same flattening satisfies both discretized meridian constraints.'^ 

If we compute matrices A, B using one meridian path and A, B using the other 
quad type, flattening, and edge the same — the change (3.14) implies that 



keeping 



(3.15) 



A = 



where Pj^^ is the SL(A^, Z) matrix 



(3.16) Pj = J -h r^Ni ■ 

i.e. the identity plus an extra entry '1' in the A^*^ (meridian) row and P^ column. Since 
detPj^^ = 1, this immediately shows that det (AA^// + BA^^) = det {AA^» + BA"^), and 
so the a change in the meridian path cannot affect r^. 

3.5. Invariance under choice of flattening. Now suppose that we choose two flattenings 
(/, /', /") and (7, f , 7"), both satisfying 

(3.17) Af + Bf" = u, / + /' + /" = !. 



(m) 



B 



(3.18) 



A/ + Bf 



/ + /' + /" = ! 



We may assume that we have a quad type with B non-degenerate. Indeed, by the result of 
Section 3.2, flattening invariance in one quad type implies flattening invariance in any quad 
type. Moreover, by Lemma A. 3 of Appendix A, a quad type with non-degenerate B always 
exists. We also note that when B is invertible the matrix B~^A is symmetric (Lemma A. 2). 

The determinant in tj- is insensitive to the change of flattening. The monomial, on the 
other hand, can be manipulated as follows. Let us choose logarithms {Z, Z\ Z") of the shape 
parameters such that AZ + BZ" = iiru. Then, assuming that B is non-degenerate, we 
compute: 



zf"z"-f 
zf"z"-f 



= exp 


[z ■ if" - 


/") - Z" ■ if - /)] 






= exp 


- Z B 


-'Aif - 7) - i^7rB~' 


u — 


B-^AZ) 


= exp 


"-iTTB"^ 


^■(/-7)] 






= exp 


: - ■ 


(/ - 7) - mB-^A/ 


if 


-/)] 


= exp 


- - zvr/" • 


if-f) + ^rrf■ir- 


J": 


] 


= exp 


\^<f" ■ 7 


-/■r)] = ±1- 







/)] 



Therefore, the monomial can change at most by a sign, and tj- is invariant as desired. 
This completes the proof of Theorem 1.3. 



□ 



Note that this would not be the case if we allowed self-intersections of the meridian loops. 
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3.6. Invariance under 2—3 moves. We finally come to the proof of Theorem 1.4, i.e.the 
invariance of tj- under 2-3 moves. We set up the problem as in Figure 7. Namely, we suppose 
that M has two different regular triangulations T and T, with and + 1 tetrahedra, 
respectively, which are related by a local 2-3 move. Let us denote the respective (triples of) 
shape parameters as 

(3.19) Z := (Xi, X2, Z3, Z^) , Z := {W,, W^, W3, Z,, Z^) . 

We fix a quad type, labeling the five tetrahedra involved in the 2-3 move as in the figure. 
We will also assume that when calculating Neumann- Zagier matrices A and B, we choose 
to ignore an edge that is not the central one of the 2-3 bipyramid. 




Figure 7. The geometry of the 2-3 move: a bipyramid split into two tetra- 
hedra for triangulation [jf^i Aj, and three for triangulation IJilt^ ^i- 

There are nine linear relations among the shapes of the tetrahedra involved in the move; 
three come from adding dihedral angles on the equatorial edges of the bipyramid: 

(3.20) W[ = Xi + X2 , W^=X[+ X'^ , = X'l + X^ , 

and six from the longitudinal edges: 

Xi = W2 + W^, X[ = W^ + W^, X'l = Wi + W!^, 

X2 = W!l + W^, X'^ = W'^ + W2, X'^ = Wl^ + Wi. 

Moreover, due to the central edge of the bipyramid, there is an extra gluing constraint in T: 

(3.22) W[ + w:^ + W^ = 2tii . 

After exponentiating the relations (3.20-3.22), and also using Ziz[z'l = —1 and z'/ + z^^ — 
1 = for every tetrahedron Aj and Aj, we find a birational map between the shape param- 
eters in the two triangulations. Explicitly, 

(3.23) 

r / / 1 - 2;^^ / 1 - a^r^ r 1 - w'2~^ 1 - ^ 

<^ = X1X2, W2 = , = \ or <xi = — —, X2 = — — \ . 

I 1 — Xi 1 — X2 J 1^ 1 " ""^3 ^ — "^2 ^ 

Note that the birational map is well defined and one-to-one as long as no shape parameters 
{xi,X2,wi,W2,W3) equal 0, 1, or 00. This condition is satisfied so long as triangulations T 
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and T are both regular. (A necessary condition is that no univalent edges are created on 
one side or the other of the 2-3 move; this is also sufficient when considering the discrete 
faithful representation of M.) 

We must also choose a flattening in the two triangulations. Let us suppose that for 
[jfj^^ Aj we have a flattening with (triples of) integer parameters / = (di, d2, d^, /s, /at). 
This automatically determines a flattening / = (ei, 62, /s, /at) for the Uj=i ^« triangula- 
tion, by simply setting 

ei = d2 + d'l , e'^ = d^ + d'[ , e'[ = di + d'^ , 

(3.24) 

62 = 4' + 4 , e'2 = d'i + d2 , e'2' = 4' + di . 

This is a discretized version of the six longitudinal relations (3.21). One can check that 
expected relations such as ci + e'^ + e'{ = 1 are satisfied by virtue of the discretized edge 
constraint d[ + 4 + = 2 (c/. (3.22)). 

We have all the data needed to calculate t-j-. Let us start with determinants. In the 
triangulation IJ^^ Aj, we write the matrices A and B schematically in columns as 

(3.25) A = (ai, a2, a^) , B = (61, 62, h) , 

with Oj meaning (as, 04, a^) and similarly for 6j. This leads to a determinant 



(3.26) 



det (AA,. + BA7I) = det (aix'/ + ^, a2x'2' + ^, a^zf + |) 



Alternatively, in the triangulation IJilt^ matrices A and B have one extra row and 

one extra column. The extra gluing condition (3.22) causes the extra row in both A and B 
to contain three — I's. Altogether, the matrices take the form 



(3.27) 
so that 



A 



bi + 62 ai «2 



B 



a2 + bi ai + 62 



AA5// + BAr^ 



—w 



-W'o - -4 



-W 



1 



(61 + 62)^1 aiW'2 



2 

a2+bi 



3 

„ „,,' I ai+b2 
^2^3 + 



{bi + b2)w[ a,w'2 + ^ 



a2w'^ + 



ai+fc2 




It is then straightforward to check, using the map (3.23), that 
(3.28) 



(AA," + BAji) 



/I 






1 








1 








V 



/I 




Vo 





1 









x'l 

1 








1/ 



AA. 



The determinant of the last matrix on the left hand side is 1 
Therefore, 



W2W^ 



BA: 



Wi. 



(3.29) 



det (AA^// + BAJ^) = -wf^ det {AA,n + BA;^) . 
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We should also consider the monomial correction. However, with flattenings related as in 
(3.24), and with shapes related by the exponentiated version of (3.21), it is easy to check 
that 

(3.30) z^"^'-^ = w,{-lf^-^"zf"z"~^ . 

We have thus arrived at the desired result; by combining (3.29) and (3.30), we find 

(3.31) det {AA,„ + BA^i) P"z"'^' = ± det (AA,» + BA~^)zf"z"'f , 

so that Tj- is invariant under the 2-3 move. This completes the proof of Theorem 1.4. □ 

4. Torsion on the character variety 

Having given a putative formula for the non-abelian torsion of a cusped hyperbolic man- 
ifold M at the discrete faithful representation po, it is natural to ask whether the formula 
generalizes to other settings. In this section, we extend the torsion formula to general rep- 
resentations p : vri(M) — > (P)SL(2, C) for manifolds M with torus boundary, essentially by 
letting the shapes z be functions of p. We also find that some special results hold when M 
is hyperbolic and the representations lie on the geometric component of the SL{2, C) 

character variety. 

We will begin with a short review of what it means for a combinatorial ideal triangulation 
to be regular with respect to a general representation p. We will also finally prove Proposition 
1.7 from the Introduction. Recall that Proposition 1.7 identified a canonical connected subset 
Xep of the set of regular triangulations Xp^ of a hyperbolic 3- manifold M. This result allowed 
us to construct the topological invariant tm- 

We then proceed to define an extended Neumann- Zagier datum (3j- = {z, A, B, /) suitable 
for a general representation p, and propose a generalization of the torsion formula: 

(4.1) Trip) ■= ±^ det (AA,,, + BA:')zf"z"'^ . 

This formula looks identical to (1.8). However, the shape parameters here are promoted to 
functions z — z{p) of the representation p, which satisfy a well known deformed version 
of the gluing equations. Moreover, the fiattening / is slightly more restricted than it was 
previously. We will prove in Section 4.5 that 

Theorem 4.1. The formula for Tf{p) is independent of the choice of extended Neumann- 
Zagier datum, and is invariant under 2-3 moves connecting p-regular triangulations. 

When M is hyperbolic, it turns out that po-regular triangulations are p-regular for all but 
finitely many representations p G Then we can create a topological invariant tm that 

is a function on just as in Proposition 1.7, by evaluating Tq-{p) on any triangulation 

in the canonical subset Xep C Xp,^. 

In general, there is a rational map from the character variety Xm to the zero-locus Ym 
of the A-polynomial AMi^i^n) = [CCG^94], for any M with torus boundary. Therefore, 
the shapes z and the torsion tj- are algebraic functions on components of the A-polynomial 
curve Ym- When M is hyperbolic and p G somewhat more is true: the shapes are 

rational function on the geometric component 5^//°™ (Proposition B.l). Then 

(4.2) TM G CiY^n = Q{mm/{Afr{^,m)) , 
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where y4|^°™(£, m) is the geometric factor of the A-polynomial. We will give a simple example 
of the function tm for the figure-eight knot in Section 4.6. 

4.1. A review of p-regular ideal triangulations. In this section we discuss the p-regular 
ideal triangulations that are needed to generalize our torsion invariant. Let M denote a 3- 
manifold with nonempty boundary and p : 7ri(M) — y PSL(2, C) a PSL(2, C) representation 
of its fundamental group. Let X denote the set of combinatorial ideal triangulations T of 
M. Matveev and Piergallini independently showed that every two elements of X with at 
least two ideal tetrahedra are connected by a sequence of 2-3 moves (and their inverses) 
[Mat87, Pie88]. For a detailed exposition, see [Mat07, BP97]. 

Given an ideal triangulation T, let Vj- denote the affine variety of non- degenerate solutions 
of the gluing equations of T corresponding to its edges. There is a developing map 

(4.3) Vt ^ Xm , 

where Xm := Hom(7ri(M), PSL(2, C))/PSL(2, C) denotes the affine variety of all PSL(2,C) 
representations of 7ri(M). 

Definition 4.2. Fix a PSL(2, C) -represent at ion of M. We say that T E X is p-regular if p 
is in the image of the developing map (4.3). 

Let Xp (Z X denote the set of all p-regular ideal triangulations of M . When M is hyper- 
bolic, let pq denote its discrete faithful representation po and let C Xm denote the 
geometric component of its character variety [Thu77, NZ85]. We then have the following 
result. 

Lemma 4.3. (a) T G Xp^ if and only if if T has no homotopically peripheral {i.e., univalent) 
edges. 

(b) If T G A'pq, then T & Xp for all but finitely many p G 

Proof. Part (a) is well known [Cha03, BDRV, Til, DGll]. For part (b), fix T G Xp^. Observe 
that T is p-regular if the image of every edge^ of T under p does not commute with the 
image under p of the peripheral subgroup of M. This is an algebraic condition on p, and 
moreover, when p G is analytically nearby po, the condition is satisfied. It follows 

that the set of points of that satisfy the above condition is Zariski open. On the 

other hand, is an affine curve [Thu77, NZ85]. It follows that T is p-regular for all but 

finitely many p G □ 

4.2. The Epstein-Penner cell decomposition and its triangulations. Now we con- 
sider the canonical ideal cell decomposition of a hyperbolic manifold M with cusps [EP88], 
and finally prove Proposition 1.7. It is easy to see that every convex ideal polyhedron can 
be triangulated into ideal tetrahedra with non-degenerate shapes, see for instance [HRS]. 
One wishes to know that every two such triangulations are related by a sequence of 2-3 
moves. This is a combinatorial problem of convex geometry which we summarize below. For 
a detailed discussion, the reader may consult the book [DLRSIO] and references therein. 



^Note that every edge can be completed to a closed loop by adding a path on the boundary T^. The 
choice of completion does not matter for studying commutation with the peripheral subgroup. 
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Fix a convex polytope P in W^'. One can consider the set of triangulations of P. When 
d = 2, P is a polygon and it is known that every two triangulations are related by a sequence 
of flips. For general flips are replaced by geometric bistellar moves. When > 5, it is 
known that the graph of triangulations (with edges given by geometric bistellar flips) is 
not connected, and has isolated vertices. For (i = 3, it is not known whether the graph is 
connected. 

The situation is much better when one considers regular triangulations of P. In that case, 
the corresponding graph of regular triangulations is connected, an in fact it is the edge set 
of the secondary polytope of P. When d = 3 and P is convex and in general position, then 
the only geometric bistellar move is the 2-3 move where the added edge that appears in 
the move is an edge that connects two vertices of P. When d = 3 and P is not in general 
position, the same conclusion holds as long as one allows for tetrahedra that are flat, i.e., lie 
on a 2-dimensional plane. 

Returning to the Epstein-Penner ideal cell decomposition, let Xep denote the set of reg- 
ular (in the sense of polytopes and in the sense of po) ideal triangulations of the ideal cell 
decomposition. The above discussion together with the fact that no edge of the ideal cell 
decomposition is univalent, implies that Xep is a connected subset of Xp^. This concludes 
the proof of Proposition 1.7. 

4.3. Neumann-Zagier datum and the geometric component. Let M be a manifold 
with torus boundary and T a (combinatorial) ideal triangulation. The Neumann-Zagier 
datum Pf = [z, A, B) may be generalized for representations p G Xm besides the discrete 
faithful. 

To begin, choose a representation p : 7ii{M) — )• PSL(2, C), and, if desired, a lift to 
SL(2, C). Let (/i. A) be meridian and longitude cycles^ on dM, and let (m='=^,£^^) be the 
eigenvalues of p(/i) and p(A), respectively. For example, for the lift of the discrete faithful 
representation to SL(2, C), we have (m, i) = (1,-1) [Cal06]. These eigenvalues define a map 

(4.4) Xm (C*) VZ2 , 

whose image is a curve Ym, the zero-locus of the A-polynomial AM{i,m) = [CCG"'"94]. (If 
we do not lift to SL(2,C), then Ym is a corresponding "PSL(2,C)" A-polynomial.) We will 
denote the representation p as pm to emphasize its meridian eigenvalue. 

Now, given a triangulation T, and with A, B, and u defined as in Section 2.3, the gluing 
equations (2.10) can easily be deformed to account for m ^ 1. Namely, we find [NZ85] 

N 

(4.5) Yl Zj^^^z'^'^'^ = {-ly^m^^^^ . 

i=i 

The developing map (4.3) maps every solution of these equations to a representation 7ri(M) — 
PSL(2, C) with meridian eigenvalue ±m. The triangulation T is pm-regular if and only if pm 
is in the image of this map. We can similarly express the longitude eigenvalue as a product 



Recall again that these cycles are only canonically defined for knot complements. In general there is 
some freedom in choosing them, but the torsion depends in a predictable way on the choice, cf. [Yam08]. 
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of shape parameters 

N 

(4.6) l[z,^''^zp''^ = {-lf'H', 

i=i 

for some 2Cj,2Dj,2ux G Z (the extra factors of two will be explained below). Then, if T 
is a Pm-regular triangulation, the irreducible component of Ym containing is explicitly 
obtained by eliminating all shapes Zj from (4.5)-(4.6). In order to lift to SL(2, C), we take 
an appropriate square root of the longitude and meridian equations. 

In general, the shapes Zj are algebraic functions on components of the variety Ym- However, 
if M is hyperbolic and T is regular for all but finitely many representations on the geometric 
component V^f/"™, then the shapes zj become rational functions, Zj G C(Y^'^°™). We provide 
a proof of this fact in Appendix B. The field of functions C(Fjf/°™) may be identified with 
Q(m)[£]/(A^'^°™(£, m)) , and the functions Zj{i,m) can easily be obtained from equations 
(4.5)-(4.6). 

4.4. Flattening compatible with the longitude. The Neumann- Zagier datum (3-j- = 
{z,A,'B), with z a solution to (4.5) for general m, may be enhanced by a flattening. We 
will need a generalized flattening that solves a discretized version of the longitude equation 
in addition to the edge and meridian equations. 

Logarithmically, the deformation of the meridian equation (2.5) reads 

TV 

(4.7) : (Gtv+m^. + G'^+,,Z: + G^^.^Z^') = 2u ; 
and we may similarly introduce a logarithmic longitude equation 

N 

(4.8) A : (Gtv+2.^. + G'^+2,Z[ + G%^2,Z':) = 1v , 

i=\ 

where^ 

(4.9) m = e'^ , £ = -e^ 

The integers GAr+2,i, G^_,_2 j, G^_,_2 « obtained by adding and subtracting shape parame- 
ters along a simple normal curve in dM that is homotopic to the longitude. 

Just as we obtained A, B, and v from the edge and meridian equations (with or without 
deformation), we may define 

(4.10) Ci = -{GN+2,i - G'^_j,2,i) 7 A = 2 (^'iV+2,i " G'^^2,i) ! '^A = 'Y ^'N+2,i ) 

1=1 

SO that 2C-Z+2D-Z" = 2v+2Tiiv\, which exponentiates to (4.6). The remarkable symplectic 
property of A and B may be extended to C and D. Namely, there exists a completion of 

^Thc sign in ^ = —e" is not important, since the above equations are really only encoding the combinatorics 
of a triangulation. However, for hyperbolic Af , it is the appropriate sign so that dividing (4.8) by two and 
exponentiating produces the correct square root of the longitude equation in a neighborhood of the discrete 
faithful representation, c/. [Cal06]. 



THE QUANTUM CONTENT OF THE GLUING EQUATIONS 



25 



(A B) to a full symplectic matrix ( c d ) ^'^ch that the bottom rows of C and D are the 
vectors C and D [NZ85]. In particular, this means that 

(4.11) A;v-^-B;v-C= 1, 

where An,'Bn are the bottom (meridian) rows of A, B. Unfortunately, in order for the full 
symplectic matrix to be canonically normalized, extra factors of 1/2 must be included in 
(4.10). 

Now, a generalized flattening (/, /', /") is an integer solution to equations (2.11a)-(2.11b) 
and GN+2,if + G'^_,_2 if + G^_,_2 if" = 0. Equivalently, it is a solution to 

(4.12) Af + Bf" = u, C-f + D-f" = u,, 
and / + /' + /" = 1. A generalized flattening always exists [Neu92]. 

4.5. Invariance of the generalized torsion. We finally have all the required ingredients 
for the generalized torsion formula. Let M be a three-manifold with torus boundary, and pm '■ 
7ri(M) —J- (P)SL(2,C) a representation with meridian eigenvalue m. Let T be a p^-regular 
triangulation of M, which exists by Lemma 4.3 at least for a dense set of representations 
on the geometric component of the character variety. Choose an enhanced Neumann- Zagier 
datum (2;, A,B,/), with z = z{pm) satisfying the deformed gluing equations (4.5) and / 
satisfying (4.12). Then, as in (4.1), we define 

TriPrn) := ±1 det (AA,. + BA;^)z^"z"-^ . 

We can now prove Theorem 4.1. 

Repeating verbatim the arguments of Section 3, it is easy to see that tj- is independent of 
a choice of quad type, a choice of an edge of T and a choice of a meridian loop. The crucial 
observation is that the equations AZ + BZ" = inu (including the meridian equation) are 
never used in the respective proofs. Therefore, deforming the meridian equation by n 7^ 
does not affect anything. For the same reason, it is not hard to see that the formula is 
invariant under p^-regular 2-3 moves, by repeating the argument of Section 3.6. 

The only nontrivial verification required is that tj- is independent of the choice of flatten- 
ing. This does use the gluing equations in a crucial way. We check it now for m 7^ 1. 

Choose logarithms (Z, Z', Z") of the shape parameters and a logarithm u of m such that 
Z + Z' + Z" = in and 

(4.13) AZ + BZ" = 2u + iTTu , 

where u denote the A^- dimensional vector (0, 0, 0, u)^. By independence of quad type 
and Lemma A. 3, we may assume we are using a quad type with non-degenerate B. Now, 
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suppose that (/, /', /") and (/, /', /") are two different generalized ffattenings. Then: 

{Z.f"-Z"-f)-{Z-f"-Z"-f) 

= Z ■ if" - n + B-\AZ - iT^v - 2u) ■(/-/) 

= Z ■ if" - J") + Z ■ A(/ - 7) - zttB-V ■(/-/)- 2B-^u ■ if - /) 

= -zvrB-V ■(/-/)- 2B^u ■ if - f) 

= ^vr(r ■ 7- / • 7') - 2B 'u ■ if - f) , 

by manipulations similar to those of Section 3.5. 

The new term 2B^^u ■(/ — /) is now dealt with with by completing the Neumann- 
Zagier matrices (A B) to a full symplectic matrix (c d) ^ Sp(2A^, Q), whose bottom row 
agrees with (C, -D). The symplectic condition implies that AD^ — BC^ = /, or B"^ = 
B-^AD^ - C^. Then 

B 'u ■ if - 7) = B"^AD^ii ■ (/ - 7) - C'u ■ if - 7) 
= u-m~'Kif-f)-u-Cif-f) 
= -u . (C(/ - /) + D(r - h) 
= -u{C-if-f) + D.if"-f")). 

In this last equation, only the bottom row of C and D appears, due to the contraction with 
u = (0,0, ...,0,u). But this bottom row is precisely what enters the generalized flattening 
equations (4.12); since both fiattenings satisfy these equations, we must have B~^'U-(/ — /) = 
0. Therefore, upon exponentiating, we find 

(4.14) z^"z"-^ = i-iy"-^'-^-^'"z^'"z"-^' = ±z^"z"-~i , 

which demonstrates that r-j is independent of the choice of flattening. Theorem 4.1 follows. 

□ 

4.6. Example: 4i continued. We briefly demonstrate the generalized torsion formula, 
using representations on the geometric component of the character variety X^^^ for the figure- 
eight knot complement. 




Figure 8. Longitude path for the figure-eight knot complement. 
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We may consider the same triangulation as in Section 2.5. The edge and meridian equa- 
tions (2.14) are deformed to 

(4.15) z^w^z"w" = 1 , zwz" = -m^ , 

with z" = \ — z~^, w" = 1 — as usual In addition, there is a longitude equation 
that may be read off from the longitude path in Figure 8. In logarithmic form, we have 
-2Z + 2Z' = 2v, or 

(4.16) - 4:Z - 2Z" = 2v - 27ri , 
from which we identify 

(4.17) C = (-2,0), D = (-1,0), z/A = -l. 
Dividing (4.16) by two and exponentiating, we find 

(4.18) z-''z"-^ = i, 

with i = —e". This is the appropriate square root of (4.6) for lifting the geometric repre- 
sentations to SL(2, C). We can easily check it: by eliminating shape parameters from (4.15) 
and (4.18), we recover the geometric SL(2,C) A-polynomial for the figure-eight knot, 

(4.19) v4|7'"(£, m) = - {I - rr? - 2m^ - m*^ + m!'y + mH^ . 

We may also use equations (4.15)-(4.18) to express the shape parameters as functions of 
I and m. We find 

m? — m'^ m? + £ 



(4.20) z = — r— ^ , w 2 -2 • 

These are functions on the curve = {y4|°°™(£, m) = 0}. 

The flattening (2.15) does not satisfy the new longitude constraint C ■/ + £)■/" = i^a, so 
we must find one that does. The choice 

(4.21) (/„ U, /:,) = (0, 0, 1; 0, 0, 1) 

will work. Repeating the calculation of Section 3.1 with the same A and B but the new 
generalized flattening, we now obtain 



, . , 1 , ,fz" + l w" + l\ 



= ±-(z"w" -l)zw 
2^ ' 

1 — m? — 2771^ — + — 2m'^i 
4.22 = ± -—: . 

2m^(m"' — m 

This is in full agreement with the torsion found by [GM08, DGLZ09]. Note that for fixed 
m there are two choices of representation pm on the geometric component of the character 
variety; they correspond to the two solutions of A^"^{i,m) = in i. 

Remark 4.4. It is interesting to observe that the numerator of (4.22) is exactly dA^^'^ /di. 
That the numerator of the geometric torsion typically carries a factor of dA^°™/di might be 

gleaned from the structure of "A-polynomials" in [DGLZ09, OS], and will also be explored 
elsewhere. 
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5. The state integral and higher loops 

Our explicit formulas for the torsion rr, as well as higher invariants ST,n-, have been ob- 
tained from a state integral model for analytically continued SL(2, C) Chern-Simons theory. 
In this section, we will consider the entire state integral, and analyze its asymptotics in order 
to re-derive the full asymptotic expansion 

(5.1) Zr{h) = h-i exp ^Sta + 5r,i + ^5r,2 + ^^5r,3 + • • • 

and to unify the formulas of previous sections. 

The basic idea of a state integral is to cut a manifold M into canonical pieces (ideal 
tetrahedra); to assign a simple partition function to each piece (a quantum dilogarithm); and 
then to multiply these simple partition functions together and integrate out over boundary 
conditions in order to obtain the partition function of the glued manifold M. A state integral 
provides a finite- dimensional reduction of the full Feynman path integral on M. 

Currently, there are two flavors of SL{2, C) state integrals in the literature. The first, 
introduced in [HikOl, Hik07], studied in [DGLZ09], and made mathematically rigorous in 
[AK], is based on a 3-dimensional lift of the 2-dimensional quantum Teichmiiller theory in 
Kashaev's formalism [Kas98]. It uses variables associated to faces of tetrahedra. The second, 
developed in [Dim], explicitly uses shape parameters — associated to edges of tetrahedra — 
and constitutes a 3d lift of Teichmiiller theory in the Fock-Chekhov formalism [FC99]. The 
two types of state integrals should be equivalent, though this has only been demonstrated 
in isolated examples so far [SV]. 

It is the second state integral that we employ in this paper, due to its explicit dependence 
on shape parameters. Indeed, suppose that M is an oriented one-cusped hyperbolic manifold 
with a po^regular triangulation T and enhanced Neumann-Zagier datum Pj- = [z, A, B, /), 
with A/ + B/" = u. We must also assume that B is non-degenerate, which (Lemma A. 3) 
is always possible. Then we will^ show in Appendix C that the state integral of [Dim] takes 
the form 

1=1 

where iph{Z) is a non-compact quantum dilogarithm [Bar99, Fad95], the Chern-Simons par- 
tition function of a single tetrahedron. The integration variables Zi are, literally, the loga- 
rithmic shape parameters of T. 

The integration contour of (5.2) is unspecified. A complete, non-perturbative definition 
of Z-T-{h) requires a choice of contour, and the choice leading to invariance under 2-3 moves 
(etc.) may be quite subtle. However, a formal asymptotic expansion of the state integral as 
in (5.1) does not require a choice of contour. It simply requires a choice of critical point for 
the integrand. Then the asymptotic series may be developed via formal Gaussian integration 
in an infinitesimal neighborhood of the critical point. 



'''Here we multiply (C.23) (at it = 0) by an extra, canonical normalization factor (27r/fi,)'^/^, in order to 
precisely match the asymptotics of the the Kashacv invariant at the discrete faithful representation. 
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We will show in Section 5.1 that all the leading-order critical points of (5.2) are logarithmic 
solutions to the gluing equations 

(5.3) critical points < — > z^{l - z'^)^ = (-1)" , 

with z = exp(Z). In particular, the critical points are isolated. Then, choosing the discrete 
faithful solution to (5.3), we formally expand the state integral to find that 

• 5'7-,0; the evaluation of leading-order part of the integrand at the critical point, is the 
complex volume of M; 

• exp(— 2S'7-,i) is expressed as the determinant of a Hessian matrix 

7^ = _B-i A + Ari, , 

with a suitable monomial correction, and reproduces the torsion (1.8); and 

• the higher S'j-^n are obtained via a finite-dimensional Feynman calculus, and explicitly 
appear as rational functions of shape parameters. 

It follows from the formalism of [Dim], reviewed in Appendix C, that the state integral 
(5.2) is only well defined up to multiplicative prefactors of the form 

/TT^ in h \ 

(5.4) len" T 24 / ' a,b,ceZ. 

This means that we only obtain (^S-j-^.tj- = 47r^e~^'^'^'i , S'7-,2) modulo (^Z, i, jg'^), respec- 
tively; however, all the higher invariants Sj-^n>3 should be unambiguous. Moreover, in Section 
3 we saw that the ambiguity in tj- could be lifted^ to a sign ibl. Although the construction 
of the asymptotic series (5.1) appears to depend on T, we certainly expect that 

Conjecture 5.1. The invariants {S'j-^n}^_Q are independent of the choice of regular trian- 
gulation and Neumann- Zagier datum (including the choice of quad type with det B 7^ 0, 
etc.), up to the ambiguity (5.4), and thus constitute topological invariants of M. 

We now proceed to analyze the critical points and asymptotics of (5.2) in greater detail. 
In Section 5.5, we will also generalize the state integral to arbitrary representations, with 
non-unit meridian eigenvalue m = e" 7^ 1, and give an example of 5*7-^2 ('^)) S-j-^sim) as 
functions on the character variety Y^j°"^ for the figure-eight knot. 

5.1. Critical points. We begin by showing that the critical points of (5.2) are indeed 
solutions to the gluing equations. For this purpose, we need to know the quantum dilogarithm 
iphiZ). It can be defined as [Fad95] 



(5.5) = n — 



1 — g^e 



r=l 



for |g| < 1, where 

, , , J -47r^ r 2'Ki 

(5.6) g:=expAi, g := exp — - — , Z := —-Z . 

lb lb 



It may also be possible to lift the ambiguities in S'7-,0 a-nd 5*7-^2 by using ordered triangulations, as in 
[Neu04] or [Zic09]. 
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iph{Z) is the Chern-Simons wavefunction of a single tetrahedron [Dim]. The quantum dilog- 
arithm has an asymptotic expansion as /i — )■ 0, given by (c/. [DGLZ09, Eqn.3.26]) 



(5.7) 



ipniZ) ~ exp^^^ — Li2_„(e^^^ 



exp 



n=0 

^Li,(e 



ni 



720 



z(l + z)z"^ + ... 



where i?„ is the n*^ Bernoulh number, with Bi = 1/2. 

The coefficients of strictly positive powers of h {i.e. > 2) in the expansion are rational 
functions of z = e^, but the two leading asymptotics — the logarithm and dilogarithm — are 
multivalued and have branch cuts. In contrast, the function ipniZ) itself is a meromorphic 
function on C for any fixed h ^ 0. Branch cuts in its asymptotics arise when families of 
poles collide in the /I — ?■ limit. In the case of purely imaginary h with lmh> (a natural 
choice in the analytic continuation of SU(2) Chern-Simons theory), a careful analysis of this 
pole-collision process leads to branch cuts for Li2 and Lii that are different from the standard 
ones (Figure 9). We indicate the modified analytic structure of these two functions (really 
functions of Z rather than e~^) with an extra tilde. 



/WWWVWVWWN4 H — 27ri 



Im(Z) 
47r'i 



Lii,2(e" 



Re(Z) 



lm{Z) 



2-Ki 



Lii,2(e" 



-2m 



-4-Ki 



Figure 9. Rotating the standard branch cuts of Li2(e ^) and Lii(e ^) to 
produce Li2{e~^) and Lii(e~^), as functions of Z. The shaded region indicates 
where the standard logarithms of shape parameters for the discrete faithful 
representation lie. 

Now, the critical points of the integrand, at leading order^ in the h expansion, are solutions 



to 



— ( - y / ■ B-V - znZ ■ B-V + -Z ■ B-'AZ + ^^2{e 



-z7r(B-V), + (B-'AZ), - Lii(e-^'; 



^Wc treat all sub-leading terms as perturbations. The exact location of the critical point will acquire 
perturbative corrections, described in Section 5.4. 
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in other words, 

(5.8) AZ + B(-Lii(e"^)) =i7ri/. 

Since exp[— Lii(e~'^^)] = 1 — we see that every solution to (5.8) is a particular logarithmic 
lift of a solution to the actual gluing equations z-^{l — z~^)^ = (—1)''. It is a lift that precisely 
satisfies the logarithmic constraints (2.6) of Section 2, with Z" = — Lii(e~^'). 

When < Im Zi < vr, the branches of the standard logarithms and dilogarithms agree with 
those of the modified ones. In particular, given the discrete faithful solution to z'^{l—z~^)^ = 
(— 1)", taking standard logarithms immediately produces a solution to (5.8). Therefore, the 
discrete faithful representation always corresponds to a critical point of the state integral. 



5.2. Volume. By substituting a solution to (5.8) back into the h ^ (leading order) part of 
the integrand, we obtain the following formula for the complex volume of a representation: 



TT 



(5.9) St,o = - — / ■ B-^i^ - tnZ ■ + -Z ■ B'^AZ + ^ Lis 
2 2 

i 

Some manipulation involving the flattening can be used to recast this as 



mod 



TT 



(5.10) 



S- 



r,o 



;iZ-znf).{Z" + znf") + J2^i,{e~'^] 



mod 



"6 



where Z" := — Lii(e^^*). It is straightforward to verify that this formula is independent of 
the choice of quad type, choice of edge of T, choice of meridian loop, choice of flattening, 
and 2-3 moves defines a topological invariant, which agrees with the complex Chern-Simons 
invariant of M. Since the complex volume in this form has already been studied at length 
in the literature, we suppress the details here. 

At the discrete faithful representation, we can remove the "tildes" from the logarithm and 
dilogarithm. If we consider the discrete faithful solution to z^(l — z~^)^ = (—1)", and take 



standard logarithms Zi = log^j, Z-' = log(l 



(with < ImZ, ImZ" < tt), we find 



(5.11) 



5r,o = i(Vol(M) -iCS(M)) 

= _l(Z-^vr/)-(Z" + ^vrr) + 5^Li2(e-^0 (mod ^) . 



This is a version of the simple formula for the complex volume given in [Neu92]. It is known 
that the ambiguity in the volume can be lifted from 7r^/6 to 27r^ using more refined methods 
[Neu04, DZ06, GZ07, Zic09]. 



5.3. Torsion revisited. Next, we can derive our torsion formula (1.8). The torsion comes 
from the hP part in the asymptotic expansion of the state integral, which has several contri- 
butions. 
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From formal Gaussian integration around a critical point (5.8), we get a determinant 
(27r^)^/2( det where 



(5.12) = ' 



2 -1 ^ 

- y / ■ B-V - inZ ■ B-V + -Z ■ B-^AZ + ^ Li2(e-^0) 



is the Hessian matrix of the exponent (at leading order h~^). Here A^/ := diag(2;^, ...,2;^), 
with 2;^ = (1 — 2;i)~^ as usual. Multiplying the determinant is the piece of the integrand, 
evaluated at the critical point. From the part of the quadratic exponential, we get 

exp (y / ■ B-V - \z . B-V) = exp (1/ ■ (B-^AZ + Z") - \z ■ (B-^A/ + 

^exp(-lz.r + V./ 

(5.13) =[z^"^'-iy^'\ 

whereas from the quantum dilogarithm at order we find 



(5.14) exp (i Lii(e^^O) =±11 ^TfT^ = ± 



Combining the determinant (27r/l)^/^ ( det "H) , the corrections (5.13)-(5.14), and the over- 
all prefactor \j -^^{^tiK)"^ 1'^ in the integral (5.2) itself, we finally obtain 
(5.15) 



^7r3 -87r3 



det B det(-B-iA + A^O det A^nzS" z"-S \ det(AA^» + BA-i)^/"^"-/ ' 
up to multiplication by a power of i; or 

(5.16) Tm := ^T^^e-'^^^ = ±^ det(AA,» + BA;^)z^"z""^ , 

just as in (1.8). Despite the fact that the original state integral only made sense for non- 
degenerate B, the final formula for the torsion is well defined for any B. 



5.4. Feynman diagrams and higher loops. The remainder of the invariants S-^-^n can be 
obtained by continuing the saddle-point (stationary phase) expansion of the state integral 
to higher order. The calculation can be systematically organized into a set of Feynman 
rules (c/. [HKK+03, Ch.9], [BIZ80, Pol05]). The resulting formulas — summarized in the 
Introduction — are explicit algebraic functions of the exponentiated shape parameters Zi, 
and belong to the invariant trace field Em- 
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To proceed, we should first re-center the integration around a critical point. Thus, we 
replace Z ^ Z + ( and integrate over (, assuming Z to be a solution to (5.8). We expand 



1 n 



[|(j7r+|) /■B-ii/-(i7r+|){Z+C)-B-i!/+|(2+C)-B-iA(Z+C)] 



2h ' ' n\ 

k=l 1=1 



In this form, the first coefficient T^'^^Z) can be identified with an overall vacuum energy, 
while the rest of the rf'\z) are vertex factors. 

Every T^'^^Z) here is a series in h, in general starting with a term. However, vf'"' must 
vanish at leading order precisely because Z is a solution to the leading-order critical point 
equations; and we have also already extracted the leading piece of t[^^ as the Gaussian 
integration measure —^C^C- Typically, 1-vertices and 2-vertices are absent from a Feynman 
calculus. Here, however, they appear because our critical point equation and the Hessian 
(respectively) are only accurate at leading order, and incur ^-corrections. (Note that the 
1-vertices and 2-vertices are counted separately in (5.19) below.) 

The vacuum energy T^^^ contributes to every S'j-,n, n > 0. Its leading-order term is 
just the complex volume (5.10), while the hP piece contains the corrections (5.13)-(5.14) to 
the torsion. At higher order in h, we have 



(5.18) r(°)(Z) = l5o + ;^%..) + ^/-B-A/ + f;^^J]Li2_.(.ri) (mod A) . 

n=2 ■ i=l 



Each Sn, n > 2, is calculated by taking the ^ part of T^^\ and adding to it an ap- 
propriate sum of Feynman diagrams. The rules for the diagrams are derived from (5.17) as 
follows. There are vertices of all valencies k = 1,2,..., with a vertex factor given by F-'^^ 
One draws all connected diagrams (graphs) with 



(5.19) # loops + # 1-vertices + # 2-vertices < n 



(k) 

Each fc- vertex is assigned a factor F) , and each edge is assigned a propagator 

(5.20) propagator : H^^- := ^H,^^ = ^(-B-^A + A^,/)r; . 

The diagrams are then evaluated by contracting the vertex factors with propagators, and 
multiplying by a standard symmetry factor. In each diagram, one should restrict to the ^""^ 
term in its evaluation. 
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Explicitly, using the asymptotic expansion (5.7) of the quantum dilogarithm, we find that 
the vertices are 

(5.21a) 1-vertex: rf ^ = --(B^M. - E ^^Li,_„(^-i) = --(R-V), + | + . . . , 

n=l 

(5.21b) 2-vertex: if = ^^Li_„(^ri) = ^ - + + . . . 

n=l 

(5.21c) fc-vertex: rf ^ = (-1)'= ^-p^Li2_„_fe(zri) (A; > 3) 

n=0 

Note that in T^^'^ we could also write B~^z/ = B~^A/ + /". When the inequality (5.19) 
is saturated, only the leading-order [hr^ or fp) terms of the vertex factors (5.21) need be 
considered. Otherwise, subleading ^-corrections may be necessary. 

Examples of 2-loop and 3-loop Feynman diagrams were given in Figures 1-3 of the Intro- 
duction, along with the entire evaluated expression for S'7-2. 

5.5. n-loop invariants on the character variety. Just as we extended the torsion for- 
mula to general representations p G Xm in Section 4, we may now generalize the entire 
state integral. The basic result for the higher invariants Sq-^n is that their formulas remain 
completely unchanged. The shapes Zi simply become functions of the representation p, and 
satisfy deformed gluing equations (4.5)-(4.6). One must also make sure to use a generalized 
fiattening whenever it occurs, just as in Section 4. 

We note that, for a hyperbolic knot complement M = S^\K, the generalized Chern- 
Simons state integral Zm{u; K) is expected to match the asymptotic expansion of the colored 
Jones polynomials JN{K;q). Specifically, one should consider the limit 

(5.22) A^^oo, n->0, g^ = e^^ = e2" fixed, 

where m = e" is the meridian eigenvalue for a geometric representation pm in the neighbor- 
hood of the discrete faithful. This is the full Generalized Volume Conjecture of [Guk05]. 

To see how formulas for the generalized invariants Sf^n, n > 0, come about, consider the 
state integral at general meridian eigenvalue m = e". From (C.23) of Appendix C, we find 

1 

(5.23) X e^' 

where u := (0,...,0,n) and D is the block appearing in any completion of the Neumann- 
Zagier matrices (A B) to ( ^ ^ ) G Sp{2N, Q), such that the bottom row D of D appears in 
the longitude gluing equation C ■ Z + D ■ Z" = v + 2Txiux (Section 4.4). Indeed, since we are 
contracting with -u, only this bottom row of D really matters in (5.23). 
The critical points of the state integral are now given by 



2tt-DB"itt+(27ri+n)/-B- 



-Z B 



(5.24) 



AZ + BZ" = 2u + iiru , 
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with Z" := — Lii(e~^). As expected, this is the logarithmic form of the deformed gluing 
equation (4.5). Thus, all critical points correspond to representations p = pm ^ Xm- The 
multivalued nature of this equation must be carefully studied to make sure desired solutions 
actually exist. However, for example, representations on the geometric component A|,j°™ 
always exist in a neighborhood of the discrete faithful representation, if we choose u = log m 
to be close to zero (and use a regular triangulation). 

We then start expanding the state integral around a critical point, setting 

(5.25) Zr{u] h) ~ ^i-iexp \\Sro{.u) + 5ri(u) + hSr2{.m) + h^Sr-i{,m) 

In 

The leading contribution S-j-^Q^p) is given, following some standard manipulations using the 
generalized flattening, by 

(5.26) Sr,o{u) = u v{u) - ^{Z - inf) ■ {Z" + mf") + J] Li2(e-^) (mod 

i=l 

Here we write Sq-fi as a function of the logarithmic meridian eigenvalue m, though a fixed 
choice of representation p will implicitly fix the choice of longitude eigenvalue v = log(— ^) as 
well. Expression (5.26) is a holomorphic version of the complex volume of a cusped manifold 
with deformed cusp. Explicitly, 

(5.27) Sr,o{u) = «(Vo1m(m) + iCSAf{u)) - 2v ^{u) . 

This is the correct form of the complex volume to use in the Generalized Volume Conjecture, 
c/. [GM08]. 

At first subleading order, we re-derive the generalized torsion formula. The calculation is 
identical to that of Section 5.3, with the exception of the correction (5.13) coming from the 
part of the exponential. This correction now becomes 

%Tl 1 

(5.28) exp yf ■ B~^w + —f ■ B^V - -Z ■ B^u 

To simplify this correction, we must use A/ + Bf" = v and a deformed gluing equation 
KZ + BZ" = 2u + iiru. The n-dependent part of the gluing equation cancels the new 
w-dependent term in (5.28), ultimately leading to the same result 



exp 

and therefore the same torsion^° 



/.B-iu+-/-B-V--Z-B-V 

2 2 



(5.29) Tr = 47re-2^^-i = - det (AA,» + BAj^) z-f" z"-^ 



1 
2 

Finally, we can produce a generalized version of the Feynman rules of Section 5.4. We 
note, however, that the w-dependent terms in (5.23) do not contribute to either the vacuum 
energy F^^^ (at order or higher), the propagator, or the vertex factors F^-'^''. Therefore, 



"'^'^The normalization of the torsion here differs from the torsion at the discrete faithful by a factor of tt^. 
In fact, we intentionally changed the normalization of the entire state integral (5.23) by tt^. This is because 
we wanted the state integral to match the asymptotics of the colored ,Iones polynomials exactly, and the 
asymptotics happen to jump by tt^ when u ^ 0, cf. [GM08]. 
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the Feynman rules must look exactly the same. The only difference is that the critical 
point equation (5.24) requires us to use shape parameters that satisfy the generalized gluing 
equations. 

5.6. Example: 4i completed. We may demonstrate the power of the Feynman- diagram 
approach by computing the first two subleading corrections Sra and Sr,z for the figure-eight 
knot complement. 

We can use the same Neumann- Zagier datum described in Section 2.5, along with the 
generalized flattening of Section 4.6. Let us specialize to representations pm on the geometric 
component of the character variety. Then the two shapes z, w are expressed as functions on 
the A-polynomial curve, 

. m^-m~'^ + i 

5.30 z = — — , w = — — , 

1 + rw^i m'^ — m ^ 

as in (4.20). 

The 2-loop invariant is explicitly given in (1.19) of the Introduction. Evaluating this 
expression in Mathematica, we find 

w^{z + l) + w'^{{ll-Sz)z-A) + w{z-l){z{z + l2)-h) + {z-2){z-lf 
~ 12{w + z-lf ■ 

Upon using (5.30) to substitute rational functions for z and w, the answer may be most 
simply expressed as 

(5.31) 2 = t = -— (m-6 - m-^ - 2m-'' + 15 - 2m^ - + m^) , 

where we have divided by a power of the torsion as suggested in (1.14) of the Introduction. 
(We have also absorbed a constant 1/8, recalling that our formula is only well defined modulo 
Z/24. 

In a similar way, we may calculate the 3- loop invariant, finding unambiguously 

(5.32) S4„3 = %^ = — (m"*' - - 2m~' + h-2rrv'-m^ + m^) . 

''41 128 

These answers agree perfectly with the findings of [DGLZ09], and the comparison there 
to the asymptotics of the colored Jones polynomials at general u. Moreover, at the discrete 
faithful representation we obtain 

^^•^^^ ^^^'^ = 7^ = "T^ ' ^^-^ = "^=128^' 

in agreement with known asymptotics of the Kashaev invariant. 

Appendix A. Symplectic properties of A and B 

The N X N Neumann- Zagier matrices A and B form the top half of a symplectic ma- 
trix (cE) e Sp(2Ar,Q) [NZ85]. In this section we discuss some elementary properties of 
symplectic matrices. 
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Lemma A.l. The x 2N matrix (A B) is the upper half of a symplectic matrix if and 
only if AB^ is symmetric and (A B) has maximal rank N . 

Proof. It is easy to see that the rows of (A B) have zero symplectic product (with respect to 
the standard symplectic form on if and only if AB^ is symmetric. In addition they span 
a vector space of rank N if and only if (A B) has maximal rank N . The result follows. □ 

Lemma A. 2. If (A B) is the upper half of a symplectic matrix and B is non-degenerate, 
then B^^A is symmetric. 

Proof. Lemma A.l implies that AB^ is symmetric, and so is (B)~^AB^((B)~^)^. □ 

It is not true in general that B is invertible. However, after a possible change of quad 
type, we can assume that B is invertible. This is the content of the next lemma. 

Lemma A. 3. (a) Suppose (A B) is the upper half of a symplectic 2N x 2N matrix. If A 
has rank r, then r columns of A and their complementary N — r columns of B is a basis for 
the column space of (A B). 

(b) There always exists a choice of quad type for which B is non-degenerate (for any fixed 
choice of redundant edge and meridian path). 

Proof. For (a) let rank(A) = r < N. Without loss of generality, we may suppose that the 
first r columns of A are linearly independent. We want to show that, together with the last 
N — r columns of B, they form a matrix of rank A^. 

If we simultaneously multiply both A and B on the left by any nonsingular matrix U G 
GL(iV, M), both the symplectic condition and the columns are preserved. (This follows from 
the fact that ( ^ u-i,t) G Sp(2A^, M).) By allowing such a transformation, we may assume 
that A takes the block form 

(A.l) A=(^Y t) 

for some A2. Similarly, we split B into blocks of size r and N — r, 

(*-^) ^ - (I: S) ■ 

Since (A B) has full (row) rank, we see that the bottom N — r rows of B must be linearly 
independent, i.e. rank(i?3 B4) = N — r. From the symplectic condition of Lemma A.l, we 
also find that B3 + B4A2 = 0, so that rank(53 B4) < rank(i?4). This then implies that B^ 
itself must have maximal rank N — r. Therefore, the last N — r columns of B are linearly 
independent, and also independent of the first columns of A; i.e. the matrix ( ^""q ) has 
maximal rank as desired. This concludes the proof of part (a). 

For part (b) let us denote the columns of A and B as and h^. A change of quad type 
corresponding to a cyclic permutation Zi 1— )■ Z[ h> Z'- h-> Zi on the i^^ tetrahedron permutes 
the 2*^ columns of A and B as {ai,bi) {hi — Oi, — flj). Therefore, given complementary 
columns of (A B) that have full rank, we can use such permutations to move all the columns 
(up to a sign) into B. □ 
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Appendix B. Shape parameters as functions on the character variety 

In this appendix, we explain why the shape parameters are rational functions on Y^j""^, 
the geometric component of the A-polynomial curve. 

Proposition B.l. With the notation of Section 4.3, we have that Zt G C(Y^j°™) for i = 
1,...,N. 

Proof. The proof is a little technical, and presented in [DG] following Dunfield [Dun99, 
Cor. 3. 2] and the appendix to [BDRV]. For completeness, we give the details of the proof 
here. 

Consider the affine variety -R(Af, SL(2, C) = Hom(7ri, SL(2, C)) and its algebrogeometric 
quotient ^m,psl(2,c) by the conjugation action of PSL(2, C). Following Dunfield from the Ap- 
pendix to [BDRV], let R{M, SL(2, C)) denote the subvariety of _R(M, SL(2, C)) x P\C) con- 
sisting of pairs (p, z) where z is a fixed point of p(7ri((9M)). Let Xm,sl{2,c) denote the algebro- 
geometric quotient of -R(M, SL(2, C)) under the diagonal action of SL(2,C) by conjugation 
and Mobius transformations respectively. We will call elements (p, z) G R{M, SL(2, C)) aug- 
mented representations. Their images in the augmented character variety X(M, SL(2,C)) 
will be called augmented characters and will be denoted by square brackets [(p, 2;)]. Like- 
wise, replacing SL(2,C) by PSL(2, C), we can define the character variety Am,psl(2,c) and 
its augmented version X M,psh{2,c)- 

The advantage of the augmented character variety Xm,sl(2,c) is that given 7 G nildM) 
there is a regular function which sends [(p, z)] to the eigenvalue of p{'~f) corresponding to 
z, using Lemma B.2 below. In contrast, in Am,sl(2,c) only the trace -|- e~^ of p(7) is well- 
defined. Likewise, in Am,psl{2,c) (resp. Aa./,sl(2,c)) only (resp. + e~'^) is well-defined. 

From now on, we will restrict to the geometric component of the character variety Xm,sl{2,c) 
and we will fix a regular ideal triangulation T. Dunfield observes in the appendix to [BDRV] 
that the developing map 

Vj > Am,PSL(2,C) 

is well-defined, 1-1 with its image inside the geometric component. The latter is an affine 
curve [Thu77, NZ85], and each shape Zi is a rational (in fact a coordinate) function of Vf. 
It follows that each shape Zi is a rational function on AAf^psL{2,c)- There is a map 

Am,SL{2,C) > -^M,PSL(2,C) 

of degree 2, which can be used to pull Zi to a rational function on Xm,sl{2,c)- In addition, 
the map 



-^M,SL(2,C) > C* X C* 

given by {p,z) H- {M,L) is of degree 1 [Dun99, Cor. 3. 2], and its image is the curve Y^°'^. 
It follows that Zi is a rational function on Y^^""^ for i = 1, . . . , A^. □ 

Lemma B.2. Suppose A = ^ ^ SL(2, C) and c 7^ 0. Then, A is an eigenvalue of A if 

and only if 2; = (A — 2d) / (2c) is a fixed point of the corresponding Mobius transformation in 
P\C). 
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Appendix C. Deriving the state integral 

In this appendix, we explain the connection between the quantization formalism of [Dim] 
and the special state integrals (5.2) and (5.23) that led to all the formulas in the present 
paper. We will first review classical "symplectic gluing" of tetrahedra, then extend gluing 
to the quantum setting and construct the state integral. There are multiple points in the 
construction that have yet to be made mathematically rigorous, which we will try to indicate. 

C.l. Symplectic gluing. The main idea of [Dim] is that gluing of tetrahedra should be 
viewed, both classically and quantum mechanically, as a process of symplectic reduction. 

Suppose we have a one-cusped manifold M with a triangulation T = {Aj}^^. Classically, 
each tetrahedron Aj comes with a phase space 

"P^Ai = {flat SL(2, C) connections on 9Aj} 
(C.l) ^ {(Z„ Zl Z'l) e C\(27rzZ) I Z, + Z', + = zvr} , 

with (holomorphic) symplectic structure 

(C.2) UJ9A, = dZ A dZ" , 

and a Lagrangian submanifold^^ 

= {flat SL(2, C) connections that extend to Aj} 

(C.3) = {e^" + 6^^-1 = 0} C VoA, 

When gluing the tetrahedra together, we first form a product 

(C.4) £x = -Cai X ■ ■ ■ Ca^ C Px = ^9Ai X ■ ■ ■ X VsAr, ■ 

The edge constraints Xj := (G/jZj + G'j^Z'^ + G^^Zf ) — 2ni from (2.4) are functions on 
the product phase space Vx, and can be used as (holomorphic) moment maps to generate 
— 1 independent translation actions tj. Recall [NZ85] that the logarithmic meridian 
and longitude holonomies {u,v) are also functions on Vx, which Poisson-commute with all 
the edges Xj, and so are fixed under these translations. Then the phase space of M is a 
symplectic quotient, 

VdM = {flat SL(2,C) connections on dM ~ T^} ^ {{u,v) e C} 
(C.5) =Vx//{ti), 

and the A-polynomial of M (more properly, components of the A-polynomial for which the 
triangulation is regular) is the result of pulling the Lagrangian Cx through the quotient, 

(C.6) £a/= "£x//(t/)" ^ {AM(e^e")=0} C Vsm ■ 

This is quite easy to check using equations (4.5) and (4.6). 



Explicitly, VdAi is a space of flat connections on a 4-punctures sphere with parabolic holonomy at the 
four punctures; while is the subspace with trivial holonomy — hence connections that extend into the 
bulk of the tetrahedron. See, e.g., Section 2 of [DGG]. 
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C.2. Quantization. Quantum mechanically, each tetrahedron has a Hilbert space "Ha^, 
a wavefunction Z^-{Zi), and a quantum operator CgAi that annihilates the wavefunction. 
The symplectic-gluing procedure extends to the quantum setting, with appropriate quantum 
generalizations of all the above operations. Roughly, one forms a product wavefunction 

(C.7) Z^{Zi,...,Z^) = Za,®---^Za^ e =^5Ai ® ■■■®HaA^, 

and restricts the product Hilbert space using — 1 new polarizations coming from the 
edge constraints. The resulting restricted wavefunction is Z;^j{u), and it is annihilated by a 
quantized version of the A-polynomial [Guk05, Gar04]. 

To make this more precise, let M again be an oriented one-cusped manifold, and choose a 
triangulation T = {Aj}^^ (regular with respect to some desired family of representations), 
a quad type, a redundant edge, and a meridian path — just as in Section 2. 

To each tetrahedron Aj we associate a boundary Hilbert space T-LoAi- It is some extension^^ 
of L^fM) that includes the wavefunction 



(C.8) z^^iz,;h) ■.= Mz^), 

where ipniZi) is Faddeev's quantum dilogarithm (5.5) [Fad95]. We also associate to Aj an 
algebra of operators 

(C.9) ^A. = C{Z„ Z[, Z'l)l{Z, + Z[ + Z'l = ^vr + I) , 

with commutation relations 

(CIO) [4^:] = [4^r] = [^;',^^] = ^- 

Then the quantization of the Lagrangian (C.3) annihilates the wavefunction, 

(C.ll) ^A, := e^" + e-^' - 1 , £a,2a. = , 

where the operators act in the representation 

(C.12) Z, = Z,, Z'l = hdzr, or e^'Z(Z,) = e^»Z(Z,) , e^'^ Z{Z,) = Z{Z,^ K) . 

In order to glue the tetrahedra together, we start by forming the product wavefunction 
Z^iZx-, Z^^ = Zai(^i) ■ ■ ■ Z^j^{Ziy). This is an element of a product Hilbert space (C.7). 
Acting on this product Hilbert space is the product Ax of algebras (C.9), which is simply 
generated by all the Zj, Z'-, Z-', with canonical commutation relations (C.IO) (and operators 
from distinct tetrahedra always commuting). 

Now, following the notation of Sections 2.2 and 4.4, we can define A^ operators Xj G Ax, 
one for each independent edge, and one for the meridian: 

'Eti [GhZ^ + G'J[ + G'/,Zf) -2m-h / = 1, A^ - 1 , 

GN+i,iZi + G^^i^jZ- + G'j^_^i_jZ-' I = N . 

Similarly, we may define an operator 

(C.14) Pn ■■= j[GN+2,iZi + G'^^2,iZi + G'^+2,iZi) 



(C.13) Xi 



""^^This space has not been mathematically defined yet; constructions of (e.g.) [AK] might prove useful 
achieving this. 
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corresponding to the longitude. Due to the symplectic structure found in [NZ85], we know 
that we may complete the set {Xi, X^v, -P/v} to a full canonical basis of the algebra Ax- 
We do this by adding N — 1 additional operators Pi, which are linear combinations of the 
Z^s, such that 

(C.15) [Pi,X,]=5i,h, [Pi,P,] = [Xj,X,] = 0, l</,j<iV. 

The operators Xj, Pj have a simple interpretation in terms of a generalized Neumann- 
Zagier datum. Namely, if we complete (A B) and the rows C, D (of Section 4.4) to a full 
symplectic matrix ( c d)' then 



(C.16) 




Here v is precisely the vector of N integers introduced in (2.8); while up = (*, *, ux), with 
ux from (4.10). The first N — 1 entries of z/p depend on the precise completion of the canonical 
basis (or the symplectic matrix), and ultimately drop out of the gluing construction. 

C.3. Quantum reduction. Classically, in order to glue we would want to set the N — 1 edge 
constraints Xj — )■ 0, and the meridian Xn — > 2u. In Section C.l, these functions were actu- 
ally used as moment maps to perform a symplectic reduction. Now we should do the same 
thing quantum mechanically. In order to reduce the product wavefunction Zx{Zi, Z^^) 
of (C.7) to the final wavefunction Zm{u) of the glued manifold M, we must transform the 
wavefunction to a representation (or "polarization") in which the operators Xj act diago- 
nally (by multiplication). In this representation, the wavefunction depends explicitly on the 
X/. The "reduction" then simply requires fixing X/ — (0, 0, 2m). Schematically, 

(C.17) Zx(Zi,...,Z;v)"™^x(^i,...,^iv) ^ ZM(w) = ix(0,. ..,0,2m). 

The transformation from to Z^ is accomplished — formally — with the Weil represen- 
tation IZ of the affine symplectic group [Sha62, Wei64]. In particular, we need 7^(a) for the 
affine symplectic transformation a in (C.16). In Section 6 of [Dim], it was discussed in detail 
how to find TZ{a) by factoring the matrix of (C.16) into generators. Then, for example, an 
"S-type" element of the symplectic group acts via Fourier transform 

whereas a "T-type" element acts as multiplication by a quadratic exponential 
(C.19) 7^ (( ^ ? )) : f{Z) ^ f{W) = eir.^"^"^ f{W) . 

Affine shifts act either by translation or multiplication by a linear exponential. 

In the present case, there is a convenient trick that allows us to find Tl{a) without de- 
composing a into generators. We assume that the block B of the symplectic matrix is 
nondegenerate, since we know we can always choose a quad type with this property. For the 
moment, let us also suppose that the affine shifts vanish, u = up = 0. Then the Weil action 
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IS 

(C.20) 



d^Z — fx-DB-iX-2ZB-iX+Z-B-iAZ' 



■ e2n' 



ZJZ). 



In particular, it can easily be verified that this correctly intertwines an action of operators 
{Zi = Zi, Z'l = hdzj on Z^{Z) with an action of operators {Xj = Xj, Pj = hdxi) on 
Zx{X). For example, 

d^Ze^(^ °^"^-^^-^"^+^-^"^^) [AZ + BZ")Z^{Z) 



7^ ('X DB-iX-2Z B-iX+Z-B-iAZl 



{AZ - hBdz)e2h^ 



Z^{Z) 



= XZ^{X). 

Non-zero affine shifts v and vp further modify the result to 



Z^{X) 



d^'Z 



^X.{^n + | 



Up 



Z -B-^AZ 



Z^{Z) 



and then, after setting X ^ 2u = (0, 0, 2u) as in (C.17), we find 



ZmM 



d'^Z 



Vd^y (27rzn)^/2 



exp 



-(27ri + fijuxu + 
n 



— ( i2u+ (ZTT+^W 



(C.21) 



Z -B-^AZ] 



BB^ (2u + (^in + fji?j-2Z- (^2w + (in 

N 



1=1 



This is the partition function of the one-cusped manifold M, modulo a multiplicative ambi- 
guity of the form exp + ^c] for a,b,c & Z, which we will say more about in Section 

C.5. By construction, this partition function is annihilated by the quantum A-polynomial 
of M. 



C.4. Introducing a flattening. In order to obtain the state integral (5.23) appearing in 
the paper, we can introduce a generalized flattening (as in Section 4.4) and use it to simplify 
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(C.22) 



(C.21). Note that the discrete-faithful state integral (5.2) follows immediately from (5.23) 
upon setting u = (0, 0, u) — 0. 

Suppose, then, that we have integers (/, /") that satisfy 

c d) (r) = (.p^ 

for some z/p whose last entry is ux- We will assume that the entire vector z/p contains integers. 
This appears to always be the case. In particular, for a one-cusped manifold M, it appears 
that the completed symplectic matrix ( ^ E ) '^^^ always be arranged to lie in Sp{2N, Z) 
after adding a suitable combination of edge equations to the longitude equation. Then, for 
example 

-uxu + V ■ DB-^w = -(C/ + D/") ■ u + (D^A/ + D^B/") ■ B"^u = / ■ B-^w , 

where we used the symplectic identities D-^B = B^D and D"^A = / + B"^C; and 

V ■ DB-V = V ■ D(/" + B-^A/) = V ■ (D/" + C/ + B-^'^/) = / ■ B-V ^ v vp 

= f ■ B-V (mod Z) 

in a similar way. These relations allow us to write the state integral (C.21) as 
(C.23) 



1 r d^z 



Nry 1 



v/d^y (27rn)^/2 



N 



i=l 



just as in (5.23). (We drop a factor of y/i from the measure, since it can be absorbed in the 
overall normalization ambiguity.) 

C.5. Normalization and invariance. The normalization of Chern-Simons state integrals 
has always been a subtle issue. For the integral of [Dim], ambiguities in the normalization 
come from two sources: the projectivity of the Weil representation, and the incomplete 
invariance of the integral (even formally) under a change of "quad type" and a 2-3 move. 

Let us consider the Weil representation first. We will assume that all symplectic matrices 
are in Sp{2N,'Z), and that all shifts involve integers (like u and up) times in + |. This as- 
sumption (which, again, is only an observed property) allows us to improve on the estimates 
of [Dim] (c/. Eqn.(6.6) there). The Weil representation becomes a projective unitary repre- 
sentation of ISp{2N, Z) ~ Sp{2N, Z) K [(z7r+ |)Z] on L^{R^), for h pure imaginary. Our 
Hilbert space l-C^^ is very close to L^(]R^), so we may hope that the Weil representation is 
also unitary projective there. The most severe projective ambiguity arises from a violation of 
expected commutation relations between shifts and T-type transformations such as (C.19). 
This leads to projective factors of the form 

(C.24) ^^P[^(^^+2j ^ =^^P[l-^ + T + 7-'^ ' 



2h\ 2, 

With the exception of factors like this, unitarity with respect to the norm = J (±2^1/0^/^ l/(^)P 
may be used to normalize Weil transformations. For example, the factor [(27ri?i)^ det B] 



"'^''A potential route to proving this statement might be similar to the flattening analysis of [Neu92]. 
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in (C.20) follows easily from formal manipulations on the integral transformation to demon- 
strate unitarity. 

The lack of complete invariance under a change of quad type (cyclic permutation invari- 
ance) and a 2-3 move can also ruin the normalization of the state integral. The change 
of quad type was analyzed, formally, in Section 6.2.1 of [Dim]. A cyclic permutation of a 
tetrahedron is accomplished by an affine version of the element ST G Sp(2N;1j), under the 
Weil representation. The single-tetrahedron wavefunction transforms as 



(C.25) MZ) ^ 



J v27ri/l 



The last equality follows from the Fourier transform of the quantum dilogarithm [FKVOl, 
PTOl]. This shows that the tetrahedron wavefunction is invariant under permutations, up 
to a factor 



(C.26) exp 



'TT in h 
.6^^T-24'" 



aeZ. 



The analysis of the 2-3 move is slightly more involved. It was done in terms of operator 
algebra in [Dim], and then explained in terms of wavefunctions in Section 6.2 of [DGG]. The 
main idea is that a 2-3 move can be done locally during the gluing procedure, by performing 
a formal, "local" transformation on the state integral. The crucial property involved is the 
Ramanujan-like identity for the quantum dilogarithm [FKVOl, PTOl], which expresses three 
quantum dilogarithms as an integral of two; for example, 

(C.27) ~ / -£=e^(^'+2^^^-(2-^+'^)(^^+^^+^))^,(-^)V^;.(^-W^i') 

J V^nih 

which holds up to a factor that is again of the type (C.26). 

Putting together all three effects, we find that we might be able to control the overall 
normalization of the state integral up to a factor of the form 



(C.28) exp 



TT^ in 1 " 
—a + —b + —c 
Qh 4 24 J 



a,b,c ^ Z . 



Appendix D. Computer implementation and computations 

An enhanced Neumann- Zagier datum is a tuple {z, A, B, /) attached to a regular ideal 
triangulation of a cusped hyperbolic manifold M. The program SnapPy [CDW] in its python 
and sage implementation computes the gluing matrices G, G', G" of Sections 2.4 and 4.4; 
and therefore it can easily compute an enhanced Neumann- Zagier datum = {z, A, B, /). 
The shape parameters z are algebraic numbers that are computed numerically to arbitrary 
precision (eg, 10000 digits) or exactly as algebraic numbers. 

A Mathematica module of the authors computes (numerically or exactly) the n-loop in- 
variants Sf^n for n = 0, 2, 3 as well as our torsion tj- given as input the Neumann- Zagier 
datum. As an example, consider the hyperbolic knot 9i2 with volume 8.836642343 . . . and 
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the SnapPy ideal triangulation with 10 tetrahedra. Its invariant trace field Eg^^ is Q(a;) 
where x = -0.06265158 + i 1.24990458 ... is a root of 

- 8x^^ + 32x^^ - 89x^^ + 195x^^ - 353x^^ + 542a:" - 719x^° + 834x^ 
- 851x^ + 764x^ - 605x^ + 421x^ - 253x^ + 130x^ - 55^^ + 18x - 3 = 

Eq^^ is of type [1, 8] with discriminant 3 ■ 298171 ■ 5210119 ■ 156953399. Our torsion is given 
by 

rg,^ = i (15 - 7a; - 15^^ + 55a;^ - 67x^ + 81x^ - 43x^ - 112x^ + 303a;^ - 488x^ 

+606a;^° - 595x" + 464a;^2 - 289x^^ + 143x^^ - 49x^^ + 8x^^) 
= -3.133657804174628986 ■ ■ ■ + 14.061239582208047255 ...i 

The two and three-loop invariants simplify considerably when multiplied by Tq^^ and Tg^^ 
respectively and are given by 

S9i2,2 rg\, = (36263 - 194718a; + 503316a;^ - 971739x^ + 1582041x^ - 2152164a;^ 

2 ■ 3 

+2372779x^ - 2109742x^ + 1426659x® - 484152x^ - 374803x^° + 836963x" 
-859483x^^ + 621288a;^^ - 326550x^^ + 109607x^^ - 16840a;^^) 
= 398.62270435384630954 • • ■ + 948.91209325049603870 ...i 



^912,3^^2 = ^ (2320213 - 19092785a;^ + 725899533;^ - 186402605x^ + 382362100a;^ 

-661985976x^ + 982969902a;^ - 1258919324x^ + 1402544816x^ 
-1359436057x^ + 1134208276x^° - 803313515x" + 473961630x^2 
-225394732x^^ + 80872920x^^ - 19104127x^^ + 2161102x^^) 
= 71793.64335382669630 ■ ■ ■ + 204530.00105728258992 . . .i 

The norm {Ni, N2, N3) = {N{Tg^^), N{S9-^^^2'ri^J, ^iSg-^^^sTg^J) of the above algebraic num- 
bers is given by 

3 ■ 298171 ■ 5210119 • 156953399 
Ni = 



N2 
N3 



217 

173137 ■ 2497646101253660962719786587619396709848261029974343408485645575954409 

2102 . 317 

1601979456387778103376978278735985249091224621424605434099548771751970874984335599199394506045984185143 

2119 



Recall that although 5*2,912 is defined modulo an integer multiple of 1/24, 5*3,912 is defined 
without ambiguity and the numerator N3 is a prime number of 103 digits. 

For a computation of the Reidemeister torsion t^j of the discrete faithful representation 
of a cusped hyperbolic manifold M, we use a theorem of Yamaguchi [Yam08] to identify it 
with 

M — 1. U=l 

Cm at 
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where cm is the cusp shape of M and Tj[|(t) G EM\t^^] is the torsion polynomial oi M using 
the adjoint representation of SL(2, C). Using the hypertorsion package of N. Dunfield (see 
[DFJ]), we can compute as follows: 

cd Genus-Comp 

sage: import snappy, hypertorsion 

def torsion (manifold, precision=100) : 
M = snappy .Manifold(manif old) 

p = hypertorsion. hyperbolic_adjoint_torsion(M, precision) 
q = p. derivative 

rho = hypertorsion. polished_holonomy(M, precision) 
z = rho . cusp_shape() 
torsion = q(l)/z. conjugate () 
return [M.nameO, torsion] 

For the above example, we have: 

sage: torsion("9_12" ,500) 

['L105002', -3.133657804174628986... + 14.061239582208047255. . .*1] 

numerically confirming Conjecture 1.8. Further computations gives a numerical confirmation 
of Conjecture 1.8 to 1000 digits for all 59924 hyperbolic knots with at most 14 crossings. 
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